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ABSTRACT 



Context. Massive clumps associated with infrared dark clouds (IRDCs) are promising targets to study the earliest stages of high-mass 
star and star cluster formation. 

Aims. We aim to determine the degrees of CO depletion, deuterium fractionation, and ionisation in a sample of seven massive clumps 
associated with IRDCs. 

Methods. The APEX telescope was used to observe the C'''0(2- 1), H'^'CO+O-l), DCO+(3- 2), N2H+(3-2), and NjD+O- 2) transi- 
tions towards the clumps. The spectral line data are used in conjunction with the previously published and/or archival (sub)millimetre 
dust continuum observations of the sources. The data are used to derive the molecular column densities and fractional abundances 
needed in the analysis of deuterium fractionation and ionisation. 

Results. The CO molecules do not appear to be significantly depleted in the observed clumps. The DCO^/HCO* and N2D+/N2H+ 
column density ratios are about 0.0002-0.014 and 0.002-0.028, respectively. The former ratio is found to decrease as a function of gas 
kinetic temperature. A simple chemical analysis suggests that the lower limit to the ionisation degree is in the range x(e) ~ 10"^ - 10"^, 
whereas the estimated upper limits range from a few 10"* up to ~ 10""*. Lower limits to x(e) imply the cosmic-ray ionisation rate of 
H2 to lie between ~ 10"" - 10"'^ s"'. These are the first estimates of x(e) and fn, towards massive IRDCs reported so far Some 
additional molecular transitions, mostly around 216 and 231 GHz, were detected towards all sources. In particular, IRDC 18102-1800 
MMl and IRDC 18151-1208 MM2 show relatively line-rich spectra with transitions some of which could be assigned to complex 
organic molecules, though the line blending is hampering the identification. The C"*0(2 - 1) transition is frequently seen in the image 
band. 

Conclusions. The finding that CO is not depleted in the observed sources conforms to the fact that they show evidence of star forma- 
tion activity which is believed to release CO from the icy grain mantles back into the gas phase. The observed degree of deuteration 
is lower than in low-mass starless cores and protostellar envelopes. Decreasing deuteration with increasing temperature is likely to 
reflect the clump evolution. On the other hand, the association with young high-mass stars could enhance and ^(^) above the levels 
usually found in low-mass star-forming regions. On the scale probed by our observations, ambipolar diffusion cannot be a main driver 
of clump evolution unless it occurs on timescales » 10* yr 
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1. Introduction 



pieces dRathborne et al. 20071 2008; JZhang et al. 2009 

TTiEr 



Beuther & Henning 2009; Wang et al. 201 



Most, if not all, in frared dark clouds (IRDCs; |Pgrault et al. 19961 
|Egan et al. 1998 1 studied so far are fragmented into clumpy 
structures. Some of the clumps within IRDCs that have 
been studied with high-resolution interferometric obser- 
vations are found to contain still denser cores, indica- 
ting that the clumps have further fragmented into smaller 



Send offprint requests to: O. Miettinen 

* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) under programmes 081.F-9823A, 
082.F-9320A, 083.C-0230A, and 085.F-9321A. APEX is a collabora- 
tion between the MPIfR, the European Southern Observatory, and the 
Onsala Space Observatory. 

** Appendices A and B are only available in electronic form at 
|http : //www ■ aanda . org | 



Massive clumps in IRDCs are promising targets to study 
the earliest stages of high-mass star and star cluster forma- 
tion (e.g., IKauffmann & Pillai 20101 and references therein). 
(Sub)millimetre dust continuum maps of IRDCs have shown 
that the clumps within them have typical radii, masses, beam- 
averaged H2 column densities, and volume-averaged H2 number 
densities of ~ 0.1 - 0.5 pc, ~ 10^ - 1 0^ Mq, > 10^- cm^^, 
and ~ 10"^ - 10^ cm"^, respectively (e.g.. lRathborne e t al. 2006) 
IParsons e t al. 20091 |Vasyunina et al. 2009) |Rygl et al. 20101 

spectral-line obser- 
the typical gas 
10 - 20 K 



Vasyunina et al. 2009[ 
[Miettinen & Harju 2010"] r 



Molecular 

vations towards such clumps have shown 
kinetic temperature to lie in the range Tki, 



' In this work, we use the term "clump" to describe a density en- 
hancement within an IRDC with a size of about ~ 0.5 - 1 pc and a mass 
of order 10^ - 10^ Mq. The term "core" is also often used for these 
objects. 
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(|Careyetal. 1998| [Teyssier et al. 2002[ ISridharan et al. 20051 
IPillai et al. 2006b'; 'Sakai et al. 2008J h ereafter SSK08 ; 
[Zhang et al. 2011 ; Devine et al. 2010 [Ragan et aLlOTT) . 
Spectral energy distributions of clumps within IRDCs 
yield dust temperature values of Jjust ~ 10 - 50 K 
(e.g., IRathborne et al. 20101 [Henning et al. 2010| l. The stai- 
formation process within some of these sources mani- 
fests itself through embedded infrared emission (e.g., 



[Chambers et al. 20091 [Ragan et al. 2009| [Henning et al. 2010|l 
outflows (fBeuther & Sridharan 2007"^ IFallscheer et al. 20091 
ISanhueza et al. 20T0l l. and H2O and CH3OH maser emission 
(IPillai et al. 2006a; Ellingsen 2006'; [Wang et al. 2 006). Hot 
cores (Beutheret al. 200 5a, 2006: Rafliborne et al. 2007, 2008) 



and UC Hn regions ( Battersby et al. 2010| l have also been 
discovered within IRDCs. These are clear signs of high-mass 
star-formation activity. On the other hand, some of the massive 
clumps are dark at mid-infared (MIR) and do not show any 
other evidence of ongoing star formation [e.g., G024.61- 
00.33 MM2 ([Rygletal.2010|l and G028.23-00.19 GLMl 



( [Battersby et al. 2010 1]. Such sources are candidates for being 
the high-mass prestellar or precluster objects and are promising 
targets for the studies of the initial conditions of massive star 
formation. 

In this paper, we aim to further constrain the phys- 
ical and chemical characteristics of massive clumps asso- 
ciated with IRDCs. We study the degrees of CO deple- 
tion and deuterium fractionation, and the fractional ionisa- 
tion in a sample of seven clumps. So far, only a few studies 
have been carried out to investigate molecular deuteration in 
IRDCs (iPiIIai et al. 20071 IChen et al. 20101 IFontani et al. 20111 
IPillai et al. 201 11 1. The level of deuteration has the potential to 
be used as an evolutiona ry indicator (e.g., [Crapsi et al. 2005 



Emprechtinger et al. 2009 IFontani et al. 201 11 1 and as such it is 
useful to study it in IRDCs. The ionisation degree, ji:(e), and the 
cosmic -ray ionisation rate of H2, ^h,. ^"e also essential knowl- 
edge towards an understanding of the initial conditions of star 
formation (e.g.. lCaseUi et al. 19981 1. Fractional ionisation deter- 
mines the coupling between the gas and magnetic field, and can 
thus play a role in the dynamical evolution of a star-forming ob- 
ject. In addition, the ionisation fraction has an influence on the 
gas-phase chemical reactions, and thus assessing the ionisation 
structure of a source helps to better understand its overall chem- 
istry. Because estimates of x(e) and ^h, in massive clumps within 
IRDCs have not been reported so far, constraining these param- 
eters through observations seems particularly worthwhile. 

This paper is organised as follows. The observations, data- 
reduction procedures, and the source sample are described in 
Sect. 2. The direct observational results are presented in Sect. 3. 
In Sect. 4, we describe the analysis and present the results of 
the physical and chemical properties of the clumps. In Sect. 5, 
we discuss our results, and in Sect. 6, we summarise the main 
conclusions of this study. 



2. Source selection, observations, and data 
reduction 

2.1. Source selection 

For this study we selected six clumps from the previous studies 
by Beuther et al. (2002a), Sridharan et al. (2005), and Rathborne 
et al. (2006), and one clump, ISOSS J18364-0221 SMMlH pre- 



viously studied by Birkmann et al. (2006) and Hennemann et al. 
(2009). We selected sources having a high mass inside a typi- 
cal clump radius of < 0.5 pc (revised masses from the above 
reference studies are ~ 60 - 360 M0; see Sect. 4.1.3), and 
which are likely to represent different evolutionary stages of 
(high-mass) star formation: two of the clumps are clearly asso- 
ciated with Spitzer 8-/im point-sources and five appear dark in 
the Spitzer 8-yL(m images. Moreover, three of the clumps are as- 
sociated with H2O and/or CH3OH masers. As will be discussed 
later, the youngest sources in our sample are likely to be GO 15. 05 
MMl and G015.31 MM3, whereas 118102 MMl appears to be 
in the most advanced stage of evolution. In this way, we were 
aiming to investigate if there are any variations or evolutionary 
trends in the clump physical and chemical properties among the 
sample. 

To gather the source sample, one criterion was that the source 
distance is not very large (revised values are ~ 2.5 - 3.5 kpc; see 
Sect. 4.1.1 and references therein) in order to achieve a reason- 
able spatial resolution with the single-dish telescope. For exam- 
ple, 20" - 30" corresponds to ~ 0.24 - 0.51 pc at the cloud 
distances, i.e., the spatial resolution is comparable to the clump 
sizes. We selected sources for which dust continuum data at 
(sub)mm wavelengths is available in order to derive the H2 col- 
umn density needed in the derivation of molecular fractional 
abundances. Also, in order to derive reliable physical parame- 
ters of the sources, they were supposed to have known gas ki- 
netic temperature or dust temperature. The source list is given in 
Table[T] and the sources are discussed in more detail in Appendix 
B. 



2.2. APEX molecular-line observations 

In this subsection we describe the molecular-line observations 
carried out with the Atacama Pathfinder Experiment (APEX) 12- 
m telescope at Llano de Chajnantor (Chile). The observed lines, 
their spectroscopic properties, and observational parameters are 
listed in Table[2] 



2.2.1. C'^0(2- 1) 

The C'^0(2 - 1) observations towards all the sources except 
J18364 SMMl were carried out on 25 March 2008 with APEX 
during the observational campaign for the Science Verification 
(SV) of the Swedish Heterodyne Facility Instrument (SHFI or 
SHeFI: Belitsky et al. 2007[ IVassilev et al. 2008al l. The source 



^ Hereafter, we use the appreviated source names, such as J18364 
SMMl etc. 



J18364 SMMl was observed on 26 April 2 009. As frontend, we 
used APEX-1 of the SHFI ( Vassil ev et aO 008b). APEX-1 ope- 
rates in single-sideband (SSB) mode using sideband separation 
mixers, and it has a sideband rejection ratio > 10 dB. The centre 
frequency for the image band is given by Vimage - Vcen ±12 GHz, 
where Vcen is the passband centre frequency, and the + refers to 
the image band corresponding to the upper or lower sideband 
(USB or LSB), respectively. The backend was the Fast Fourier 
Transfrom Spectrometer (EFTS: , Klein et al. 200 6) with a 1 GHz 
bandwidth divided into 16 384 spectral channels during SV and 
8 192 channels during observations towards J18364 SMMl. 

The observations were performed in the wobbler-switching 
(WS) mode with a 150" azimuthal throw (symmetric offsets) and 
a chopping rate of 0.5 Hz. The telescope pointing and focusing 
were checked by continuum scans on the planet Jupiter and the 
pointing was found to be accurate to ~ 3". Calibration was done 
by the chopper-wheel technique, and the intensity scale given by 
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Table 1. Source list. 



Source 


Q'2000.0 

[himis] 


"2000.0 

[°:':"] 


7 
I 

L J 


h 
U 

[°] 


rkncl 


d 


Spitz^f 8 /im/24 /jm 


IRDC 18102-1800 MM 1 


18 13 11.0 


-17 59 59 


12.624 


-0.016 


5.8 


7 7+0.5 
^- '-0.6 


point/point 


G015.05+00.07 MMl 


18 17 50.4 


-15 53 38 


15.006 


0.009 


5.9 


2.6+"-^ 


dark/dark 


IRDC 18151-1208 MM2 


18 17 50.4 


-12 07 55 


18.319 


1.792 


5.9 


-fl-5 


dark/-" 


G015. 31-00.16 MM3 


18 18 45.3 


-15 41 58 


15.281 


-0.093 


5.6 


3.0+"'* 


dark/dark 


IRDC 18182-1433 MM2 


18 21 14.9 


-14 33 06 


16.577 


-0.081 


5.1 


-'■•'-0.4 


point/point 


IRDC 18223-1243 MM3 


18 25 08.3 


-12 45 27 


18.605 


-0.075 


5.2 


3.5+°' 

^••'-0 4 


dark/point 


ISOSS J 18364-0221 SMMl 


18 36 36.0 


-02 21 45 


29.135 


2.220 


6.3 


dark/point 



Notes. Columns (2)-(7) give the equatorial and galactic coordinates [(a, (5)2ooo.o and (/, b), respectively], the source galactocentric distance 
(Rqc), and the near kinematic distance from the Sun (rf). The distances are derived by following Reid et al. (2009) as described in Sect. 4.1.1. 
Because IRDCs are seen as dark extinction features against the bright Galactic MIR background, it is expected that these clouds lie at the near 
kinematic distance. In the last column we give the comments concerning the source appearance in the Spitzer 8 and 24-pm images. 
No 24-iim image available. 



the system is T"^, the antenna temperature corrected for atmo- 
spheric attenuation. 

The spectra were reduced using the CLASS90 programme of 
the IRAM's GILDAS software packag^ The individual spectra 
were averaged and the resulting spectra were Hanning-smoothed 
in order to increase the signal-to-noise (S/N) ratio of the data. A 
first- or third-order polynomial was applied to correct the base- 
line in the resulting spectra. In one case (G015.05 MMl), a poly- 
nomial baseline of order six had to be applied because there ap- 
peared to be some emission in the wobbler off -beams which re- 
sulted in a very poor baseline. The resulting Icr rms noise levels 
are 40-60 mK at the smoothed resolutions. The final velocity 
resolution of the smoothed spectra is 0.16 km s"', except for the 
spectrum towards J18364 SMMl, where it is 0.32 km s"'. 

The 7 = 2-1 transition of C'^O contains nine hyperfine (hf) 
components. We fitted the hf structure of the transition using 
"method hfs" of CLASS90 to derive the LSR velocity (vlsr) of 
the emission, and FWHM linewidth (Av). The hf-line fitting can 
also be used to derive the line optical thickness, t. However, in 
all spectra the hf components are blended together and thus the 
optical thickness cannot be reliably determined. For the rest fre- 
quencies of the hf components, we used the values from Ladd 
et al. (1998; Table 6 therein). The adopted central frequency, 
224714.199 MHz, is that of the Jf = 29/2 -> I7/2 hf compo- 
nent which has a relative intensity of 7?, = | . 



2.2.2. Hi^^C0+(3 - 2) and DC0+(3 - 2) 

The H'^CO+(3 - 2) and DCO+(3 - 2) observations were car- 
ried out on 9-11 April and 3 June 2010. As frontend, we used 
SHFI/APEX-1 and the backend was the FFTS with a 1 GHz 
bandwidth divided into 8 192 channels. The observations were 
performed in the WS mode with similar settings as in the case of 
C'^0(2 - 1) observations. The telescope pointing and focusing 
were checked by observing the planets Saturn and Jupiter and 
the pointing was again found to be accurate to ~ 3". 

The individual spectra were averaged, smoothed, and a first- 
or third-order polynomial was applied to correct the baseline. 
The resulting Icr rms noise values are 20-80 mK in the case of 
Hi3CO+(3 - 2) and 30^0 mK in the case of DCO+(3 - 2). The 
final velocity resolutions of the smoothed H'^CO+(3 - 2) and 
DCO+(3 - 2) spectra are 0.28 and 0.34 km s"', respectively. As 
shown in the last column of Table |2] the on-source integration 
times were very different between different sources. This is par- 



and DCO+ 



Schmid-Burgk et al. 2004 



ticularly the case for Hi^CO+(3 - 2) where fint varies by a factor 
of ~ 11. This explains the large variation (factor of 4) in the 
noise values of this transition. 

The spectral-line transitions of H^^CO 
also have hf structure (e.g., 
Caselli & Pore 20051 1 . The 7 = 3-2 transition we have 
observed is split up into eight and five hf components in the 
case of H'^C0+ and DCO+, respectively. To fit the hf struc- 
ture, we used the rest frequencies of the H'^CO+(3 - 2) and 
DCO+(3 - 2) lines from the Cologne Database for Molecular 
Spectroscpy (CDMS: iMUller et al. 200 5 Fl The adopted central 
frequency of H1^C0+(3 - 2), 260255.35250 MHz, is that of 
the 7f,_f = 37/2,4 — > 25/2,3 hf component which has a relative 
intensity of ^ 0.322. We note that this frequency is 13.5 
kHz higher than the value determined by Schmid-Burgk et 
al. (2004; Table 3 therein). In the case of DCO+(3 - 2), the 
central frequency used, 216112.57790 MHz, is that of the 
7f = 34 — » 33 hf component; its relative intensity is Ri - |. 
Caselli & Dore (2005; their Table 5) determined a 2.1 kHz 
higher frequency for DCO+(3 - 2). The frequency intervals of 
the hf components for H'^CO+(3 - 2) and DCO+(3 - 2) are very 
small. Therefore, the lines overlap significantly which causes 
the hf structure to be heavily blended. 



2.2.3. N2H+(3 - 2) and N2D+(3 - 2) 

The N2H+(3 - 2) observations towards three of our target posi- 
tions (118102 MMl, G015.05 MMl, and G015.31 MM3) were 
carried out on 29 April 2009. The N2D+(3 -2) observations were 
performed during 16-17 and 19-20 October 2008, 26 and 28 
April 2009, and 7 June 2010. As frontend, we used SHFI/APEX- 
2 and APEX-1 for N2H+(3 - 2) and N2D+(3 - 2), respectively. 
The backend was the 1 GHz FFTS with 8 192 channels. The ob- 
servations were performed in the WS mode with the similar set- 
tings as described above. The pointing and focusing were done 
by continuum scans on the planets Jupiter, Mars, and Neptune. 
The pointing was found to be accurate to ~ 3" - 5". 

The data were reduced as described above, and a first- or 
third-order polynomial baseline correction resulted in Icr rms 
noise values of 40-100 mK and 20-50 mK in the case of N2H+ 
and N2D+, respectively. The final velocity resolutions of the 
smoothed N2H+ and N2D+ spectra are 0.26 and 0.32 km s 
respectively. Again, the on-source integration times were very 



Ihttp : //www ■ iram. fr/IRAMFR/GILDAS] 



|http : //www ■ astro . uni -koeln . de/cdms/catalog | 
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different for different sources, particularly for N2D^(3 - 2) [see 
Col. (10) of Table 0. 

The 7 = 3-2 transitions of both N2H+ and N2D+ contain 38 
hf components. The hf lines were fitted using the rest frequencies 
from Pagani et al. (2009b; Tables 4 and 10 therein). The adopted 
central frequencies of N2H+(3-2) and N2D+(3-2), 27951 1.832 
and 231 321.912 MHz, ai-e those of the Jf^f = 345 -> 234 hf 
component which has a relative intensity of Ri - Also in 
these cases, the hf components are blended and thus the value of 
T cannot be reliably determined through hf fitting. 



2.3. Archival dust continuum data 

Submillimetre and/or millimetre dust continuum maps of our 
sources are published in the papers by Beuther et al. (2002a; 
MAMBO 1.2 mm), WilHams et al. (2004; SCUBA 450 & 
850 iJ.m), Rathborne et al. (2006; MAMBO-II 1.2 mm), and 
Birkmann et al. (2006; SCUBA 450 & 850 yum). For the sources 
118102 MMl, 118182 MM2, and J18364 SMMl we retrieved 
850 fxm dust continuum data from the JCMT/SCUBA archive 
dPi Francesco et al. 2008 Fl No archival SCUBA data was avai- 
lable for the rest of our sources. Thus, for the sources G015.05 
MMl, G015.31 MM3, and 118223 MM3 we extracted Bolocam 
1 . 1 -mm maps from images produced by the Bolocam Galactic 
Plane Survey (BGPS; see, e.g., [Bally et al. 2010| and references 
therein^ For 118151 MM2, no archival (sub)mmdata was found 
and thus we used the MAMBO 1 .2-mm map of this source pub- 
lished in the paper by Beuther et al. (2002a; kindly provided by 
H. Beuther). 

In the present paper, the (sub)mm dust continuum data are 
used in particular to determine the H2 column densities towards 
the line-observation positions. These are used to derive the frac- 
tional abundances of the observed molecules. 



2.4. Spitzer infrared data 

In this study we also use Spitzer IRAC 8-/im and MIPS 
24-yum images of the sources to investigate the clump as- 
sociations with protostellar activity. The 8-;um images are 
from the Spitzer Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire (GLIMPSE; Benjami n et al. 200 3 Q The 24-iJ.m 
images (not shown in the paper) were taken from the Spitzer 
MIPS GALactic plane survey (MIPSGAL; [Carey et al. 2 009B 
In Fig.[T] we show the Spitzer/IRAC 8-jum images of the clumps 
overlaid with contours of (sub)mm dust continuum emission. We 
note that better quality Spitzer data of the sources are available, 
and already published (see references in Appendix B). However, 
in the present work we do not use Spitzer data to derive the phy- 
sical properties of the sources; the data are used to illustrate the 
source appearance at MIR wavelengths. 



3. Observational results 

3.1. Spectra 

The Hanning-smoothed spectra are shown in Fig. [2| Note that 
all data are presented in units of T"^. The C'^0(2 - 1) line is 
clearly detected towards all sources but the hf structure of the 
line is completely blended in all cases except GO 15. 31 MM3 



^ [http://cadcwww.dao.nrc.ca/jcmt/ 
^ fhttp : //irsa . ipac . caltech . edu/data/BOLOCAM_ GPS/' 
' http : //irsa . ipac . caltech . edu/data/SPITZER/GLIMPSE/ 
^ Ihttp : //irsa . ipac . caltech . edu/data /SPITZER/MIPSGAL/ . 



where it is partially resolved. The non-resolution of the hf struc- 
ture is caused by the linewidths to be larger than the separation 
in velocity of individual hf components. 

The H'^C0^(3 - 2) line is also clearly detected towards all 
target positions but towards GO 15. 31 MM3 the line is very weak 
(only ~ 3cr). The "double peak"-like profiles in some of the 
smoothed H'-'CO^(3 - 2) spectra can be attributed to the noise 
in most cases. In the case of J18364 SMMl, the self-absorption 
dip is a distinctive feature also in the unsmoothed spectrum, and 
thus likely to be real. DCO^(3 - 2) emission is seen towards all 
sources, but the line is very weak in 118102 MMl and GO 15. 31 
MM3 (detected only at ~ 3cr). 

N2H^(3 - 2) line-observations were performed only towards 
three of our sources, and the line is quite strong in all cases. 
The N2H+(3 - 2) line of 118102 MMl exhibits red asymmetry 
with a central dip near the velocity of the strongest hf compo- 
nent; this feature is also clearly visible in the unsmoothed spec- 
trum. The N2D^(3 - 2) line is detected towards all sources ex- 
cept G015.31 MM3, 118182 MM2, and J18364 SMMl in which 
cases the spectra have a low S/N ratio. 

3.1 .1 . Other line detections 

All the observed sources show additional spectral lines in the 
observed frequency bands, mostly in the spectra covering the 
DCO+(3 - 2) and N2D+(3 - 2) Unes. The Une identification 
was done by using Weeds, which is an extension of CLASS 
jMaret et al. 201 11 1, and applying the Jet Propulsion Laborator}0 
(JPL; IPickett et aT 1998) and CDMS spectroscopic databases. 
We used the LTE modelling application of Weeds to check if 
all predicted lines of a candidate molecule are present in the ob- 
served spectrum (see Sect. 4.2). In some cases, we were able to 
reject some line candidates on the basis of non-detection of other 
transitions expected at nearby frequencies. 

118102 MMl and 118151 MM2 show relatively Hne-rich 
specti-a at ~ 231 GHz and ~ 216 GHz, in the USB and LSB, 
respectively. Two additional lines at ~ 231 GHz (USB) are also 
visible in the spectrum of 118182 MM2 (see Fig.O. The only 
unambiguous line identification in the spectrum of 1 1 8 1 02 MM 1 
is '''CS(5 - 4). The second strongest line (T^ - 0.12 K) can be 
assigned to C'^0(2 - 1) at ~ 219.56 GHz seen in the image side 
band (marked with "image"). The strongest line in the spectrum 
(r; = 0.13 K) could be due to CNCnO{JK„,K,^ = 23o,23 - 22o,22) 
at ~ 219.53 GHz arising from the image band, or it could rep- 
resent an additional C'^0(2 - 1) velocity component along the 
line of sight. The latter is likely to be true due to the relative 
strong line intensity. The line at ~ 231 .06 GHz can be assigned 
to either OCS(19 - 18) or CH3NH2-E(7^^,;f^, = 72,5 - 7i,5) be- 
cause these two transitions are blended. However, it could also 
be caused by yet another velocity component of C'**0(2 - 1) 
seen towards 118102 MMl. From the spectrum towards 118151 
MM2 we identified the C2D(A^ = 3-2) line separated into 
two hf groups. This transition is split up into 13 hf compo- 
nents of which the detected two groups contain 10 compo- 
nents together. The hf group at ~ 216.43 GHz is blended with 
Cn2CHC''NiJK,„K^ = 252,24 - 25o,25). The Jk,„k,. = 33,o - 22,1 
transition of orf/!o-c-C3H2 is seen towards 118151 MM2, but 
the line at ~ 216.45 GHz is unidentified (marked with "U"). 
The CH3COCH3-EA(y;f„,,r^ = 159,7 - 148,6) line in the USB is 
possibly detected towards 118182 MM2. However, it is blended 
with C"*0(2 - 1) entering via the image band. The relatively 
strong intensity of = 0.17 K suggests the line is due to 
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Table 2. Observed spectral-line transitions and observational parameters. 



Transition 


V 






HPBW 


'7MB 


T 

-* sys 


Channel 


spacing" 






[MHz] 


[K] 




["] 




[K] 


[kHz] 


[kms-'] 


[min] 


DCO+(3 - 2) 


216 112.57790* 


20.7 


1.8 


28.9 


0.75 


150-184 


122.07 


0.17 


5.5-13.5 


C"0(2- 1) 


224 714.199" 


16.2 


0.01 


27.8 


0.75 


163-218 


61.04'' 


0.08'' 


5-7.5 


N2D+(3 - 2) 


231 321.912'' 


22.2 


1.7 


27.0 


0.75 


172-239 


122.07 


0.16 


6.5-48 


H'3C0+(3 - 2) 


260 255.35250* 


25.0 


3.1 


24.0 


0.74 


232-243 


122.07 


0.14 


2.5-27 


N2H+(3 - 2) 


279 511.832'' 


26.8 


3.0 


22.3 


0.74 


165-175 


122.07 


0.13 


1.5-6.5 



Notes. Columns (2)-(4) give the rest frequencies of the observed transitions (v), their upper state energies (E^/k^, where kg is the Boltzmann 
constant), and critical densities. Critical densities were calculated at T ~ 15 K (typical value in IRDCs) using the collisional rate data available 
in the Leiden Atomic and Molecular Database (LAMDA; http : //www. strw. leidenuniv.nl/~moldata/ 1 (Schoi er et al. 20051 . ForN2D+, 
we used the Einstein A-coefhcient from Pagani et al. (2009b) and the same collisional rate as for N2H^. Columns (5)-(10) give the APEX 
beamsize (HPB'W) and the main beam efficiency (77mb) at the observed frequencies, and the SSB system temperatures during the observations 
(rsy.s in scale, see text), channel widths (both in kHz and km s"') of the original data, and the on-source integration times per position (Jim). 
*'"' The original channel spacings. The final spectra were Hanning smoothed which divides the number of channels by two. From the CDMS 
spectroscopic database (Miiller et al. 2005 ). <''' From Ladd et al. (1998). *'" For the C''0(2 - 1) observations towards J18364 SMMl, channel 
width is 122.07 kHz or 0.16 km s"'. From Pagani et al. (2009b). 



C"*0(2 - 1). The weaker nearby line could be related to an ad- 
ditional C"*0(2 - 1) velocity component. 

The Jk„.k,. - 33,0-22,1 transition of orf/!o-c-C3H2 at ~ 216.3 
GHz in the" LSB is also detected towards 118102 MMl, 118182 
MM2, and 118223 MM3 (Fig.|4]i. Thermal (v = 0) SiO(6 - 5) 
line at ~ 260.5 GHz in the USB is seen towards 118102 MMl, 
118151 MM2, and 118223 MM3 (Fig. E). In particular, the 
SiO(6 - 5) line of 118151 MM2 exhibits very broad wing emis- 
sion. Finally, as in the cases of 118102 MMl and 118182 MM2, 
we detected the 7 = 2 - 1 transition of C"^0 at ~ 219.56 GHz 
in the image side band towards G015.05 MMl, 118151 MM2, 
G015.31 MM3, 118223 MM3, and J18364 SMMl (Fig.|6]l. The 
intensities of these C'**0 lines, ~ 0.09 - 0.31 K, are mostly 
comparable to those possibly detected towards 118102 MMl and 
118182 MM2. Note that the lines "leaking" from the rejected im- 
age sideband are heavily attenuated by the sideband filter. Thus, 
we cannot establish the coiTect intensity scale for the image- 
band lines. 

In Table [3] we list the observed extra transitions, their 
rest frequencies (observed frequency for the U-line), and 
upper-state energies. The rest frequencies and upper-state en- 
ergies were assigned using JPL, CDMS, and Splatalogu43 
( Reimjan et al. 2007) l spectroscopic databases. The additional 
lines/species are discussed further in Appendix A. 

3.2. Spectral line parameters 

The spectral line parameters are given in Table H) In this table 
we give the radial velocity (vlsr), FWHM linewidth (Av), peak 
intensity (T^), integrated intensity {J T^dv), peak optical depth 
of the line (tq), and excitation temperature (Tex). The values of 
Vlsr and Av for C'^O, Hi3C0+, DCO+, N2H+, N2D+, and C2D 
were derived through fitting the hf structure. For the other lines 
these parameters were derived by fitting a single Gaussian to the 
line profile. The values of were also determined by fitting 
the lines with a single Gaussian. The velocity-integrated inten- 
sities J r^dv were calculated over the velocity range given in 
Col. (6) of Table |4] The uncertainties reported in Vlsr, Av, T^, 
and J T^dv represent the formal \cr errors determined by the 
fitting routine. We used RADEf3 dvan der Tak et al. 2007] l to de- 
termine the values of tq and Tex for the lines of C^^O, H'-'CO^, 



2.312 10 2.3115 10 2.311 10 2.3105 10 



http://www.splatalogue.net7] 
" |http : //www . strw . leidenuni v . nl/ ~moldata/radex . html | 



IRDC 18102-1800 MMl 



CNCH0(23„,„-22„,„) (image) ^i.^^^-l) (image; 



- C'°0(2-1) (image) 



"cS(5-4) 
I I 



CHjNH2(72,5-7, ,,) - 
or 0CS(19-18) - 
or C"0(2-1) (i) ? 



Velocity (km/s) 
2,164 10^ 2,163b 10^ 



IRDC 18151-1208 MM! 
C2D(35/2,.-23/2,.) or 
CH2CHc'\(252^24-25o^25) ? 


^20(37/2.. -25/2,,) ^_ 

.iiliikillifiiiiiil 




-2,.,)' 


"216452.51 1 Mh| j 1 ^ ^ 




II 




Fig. 3. Multi-line spectra of 118102 MMl, 118151 MM2, and 
118182 MM2. On the spectrum towards 118151 MM2, the green 
line shows the hf-structure fit to C2D(3 - 2), and the red line 
indicates single Gaussian fits to the other lines. The relative ve- 
locities of the C2D(3 - 2) hf components are labelled with a 
short bar. On the other spectra, the green line represents single 
Gaussian fits. 



DCO+, N2H+, C-C3H2, OCS, '^CS, and SiO. RADEX modelling 
is described in more detail in Sect. 4.2 where we discuss the de- 
termination of molecular column densities. For N2D^(3 - 2) we 
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Table 3. Other candidate species/transitions observed. 



Species/ 


V 




transition" 


[MHz] 


[K] 


ortho-c-C^H2(JKa,Kc = h,Q - 22,1) V = 


216 278.7560 (JPL) 


19.5 


^I^^V^IF - ^7/2,9/2 - ^5/2,7/2i 


216 3 /2.8jUU (LUMS) 


20.8 


CH.CHC'NiJ ii„K, = 252.24 - 25„.25r 


216428.6683 (JPL) 


151. 8 


U 


216452.5110 




CNCHO(Jk„Kc = 23o.23 - 220,22)'' 


219531.6000 (CDMS) 


127.2 


c"*0(y = 2-i)'' 


219 560.3600 (JPL) 


15.8 


CH3NH2-E(yx„,x^ =72.5-7,.5)' 


231060.6041 (JPL) 


75.6 


OCS(y = 19 - 18) 


231060.9830 (JPL) 


110.9 


CH3COCH3-EA(yx„,x,, = 159.7 - 14s,6) v = 


231 080.9647 (JPL) 


95.9 


"CS(y = 5-4) 


231220.9960 (JPL) 


33.3 


SiO(y = 6-5) 


260 518.0200 (JPL) 


43.8 



Notes. For asymmetric top molecules, Ka and Kc refer to the pro- 
jection of the angular momentum along the a and c principal axes. The 
CHzCHC'^N and CNCHO transitions ai'e of type a (AK^ = 0, ±2, . . . 
and A/Cj. = ±1,±3, . . .), the 0-C-C3H2 and CH3COCH3-EA transitions 
are Z;-type (AKa = AKc = ±1,±3, ...), and the CH3NH2-E transition 
exbits c-type selection rules (AA'^ = ±1, ±3, . . . and AK^ = 0, ±2, . . .) 
l |Gordy & Cook 1984| l. This is the strongest hf component and it is 
blended with the J, F = 7/2,5/2 - 5/2,3/2 and J,F = 7/2,7/2 - 
5/2,5/2 hf components. Blended with the C2D(3 - 2) hf group. 
<''> Seen in the image band. The candidate lines of CNCHO, CH3NH2- 
E, OCS, and CH3COCH3-EA could be blended with C'*0(2 - 1) com- 
ing from the image band. The corresponding vibrational level is E] , 
/ = -1. This line is blended with 0CS(19 - 18). 



used as Tex the value derived for N2H^(3 - 2) or, when not pos- 
sible, we adopted rex[N2H+(l - 0)] from SSK08. We note that 
the Tex values obtained for the 7 = 3-2 transition of NiH^ are 
comparable to those of 7 = 1-0 from SSK08. 



4. Analysis and results 

4.1. Revision of clump properties 

In this subsection, we present recalculations of several clump 
properties presented previously in the literature (see Tables[T]and 
0. 

4.1.1. Kinematic distances 

SSK08 calculated the near kinematic distance for all of our 
clumps except J18364 SMMl (see their Table 4). They used the 
LSR velocity of N2H^(1 - 0) and the rotation curve of Clemens 
(1985) with the standard rotation parameters (0o, /?o)=(220 
km s ', 8.5 kpc), where ©o is the circular orbital speed of the 
Sun, and Rq is the galactocentric distance of the Sun. Because 
N2H^ is a tracer of high density gas, the radial velocity de- 
rived from N2H^ lines is suitable to determine the source kine- 
matic distance. The hf structure of N2H^(1 - 0) was also re- 
solved (to some degree) towards many of the sources by SSK08, 
and thus the hf-structure fitting is expected to yield a reliable 
centroid velocity. Towards the three sources for which we de- 
tected N2H^(3 - 2) lines, the obtained values of Vlsr are com- 
parable to those from N2H^(1 - 0). We recalculated the dis- 
tances from SSK08 using the recent rotation curve of Reid et al. 
(2009) which is based on measurements of trigonometric paral- 
laxes and proper motions of masers in high-mass star-forming 
regions. The best-fit rotation parameters of Reid et al. (2009) 
are (©o, /?o)=(254 km s"', 8.4 kpc). The resulting galactocentric 
distances and near kinematic distances {Rcc and d) are given in 
Cols. (6) and (7) of Table[T] The revised distances differ at most 
by 0.2 kpc from those reported by SSK08. 



The source J 18364 SMMl is an exception. Its CO(l - 0) ra- 
dial velocity of ~ 33 km s"' corresponds to a kinematic distance 
of ~ 2.2 kpc according to the Brand & Blitz (1993) rotation 
curve (see Bi rkmann et al. 2006] l. We recomputed this distance 
from the LSR velocity of C'^0(2 - 1) (34.8 km s^') because it 
is tracing higher density gas than the main CO isotopologue (si- 
milar radial velocity was obtained for the other detected transi- 
tions). The revised distance is ~ 2.5 kpc according to the Reid et 
al. (2009) rotation curve; the previous value of 2.2 kpc is within 
the errors. 

We note that for the sources 118102 MMl, 118151 MM2, 
118182 MM2, and 118223 MM3 there are kinematic distance es- 
timates (2.6, 3.0, 4.5, and 3.7 kpc, respectively) based on the 
CS(2 - 1) velocity and the rotation curve of Brand & Blitz 
(1993) (see lSridharan et al. 20 02). For both G015.05 MMl and 
G015.31 MM3, Rathborne et al. (2006) used the kinematic dis- 
tance of 3.2 kpc derived from the '^CO(l - 0) velocity and the 
rotation curve of Clemens (1985). Our distances differ by 0.1- 
1.0 kpc from the above values. 



4.1 .2. Gas and dust temperatures 

Sridharan et al. (2005) determined the NH3 rotational tempera- 
ture, r,.ot(NH3), for three of our sources (118151 MM2, 118182 
MM2, and 118223 MM3) at the 40" angular resolution. SSK08 
have determined riot(NH3) for all of our sources except J 18364 
SMMl (see their Table 10). The values of rrot(NH3) from the 
above two studies are otherwise similar except in the case of 
118223 MM3. Sridharan et al. reported the value 32.7 K for this 
source, whereas SSK08 obtained a factor of two lower tempe- 
rature (16.2^g g K) at about 1.8 times poorer angular resolution 
(73"). The ri-ot(NH3) value for J18364 SMMl at the 40" reso- 
lution is 10.75 K (Krause 2003). We converted ri-ot(NH3) from 
Krause (2003) and SSK08 into the gas kinetic temperature, r^n, 
using the relationship given by Tafalla et al. (2004; Appendix B 
therein). The uncertainty in r^in was calculated by propagating 
the larger of the two Tiot-errors given by SSK08. 

For three of our sources, G015.05 MMl, 118223 MM3, and 
J18364 SMMl, dust temperature estimates have been made by 
Rathborne et al. (2010), Beuther et al. (2010), and Birkmann 
et al. (2006), respectively. The determined dust temperature va- 
lues of G015.05 MMl, 11.0-36.0 K, bracket the r^n value of 
17.2 K. Dust temperature of 118223 MM3 was also derived by 
Beuther & Steinacker (2007) from the source SED using the 
Spitzer/MlPS data at 24 and 70 jum, and MAMBO 1.2 mm and 
PdBI 3.2 mm data. The value 15 K they obtained for the cold 
part of the spectrum is three kelvins lower than the value 18 K 
derived recently by Beuther et al. (2010). The latter also used 
the He rschel fPACS (70, 100, and 160 fim) and SPIRE (250, 350, 
and 500 fim) data, and SCUBA 850 yum data to construct the 
source SED (but not the 3.2-mm flux density). The dust tempe- 
rature of 118223 MM3 is comparable to its gas temperature of 
18.7 K. By utilising the far-infrared and submm flux density ra- 
tios, Birkmann et al. (2006) deduced the value Tdust = 16.5;^3q 
K for J 1 8364 SMMl, which is higher than the gas temperature 
11.4 K. 



The gas and dust temperatures of the clumps are given in 
Cols. (2) and (3) of Tabled 
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Table 4. Spectral line parameters. 



Source 


Transition 




Av 






/, 

To 


Tex'' 






[km s"'] 


[km .s"'] 


[K] 


[Kkms-'] 




[K] 


IRDC 18102-1800 MM 1 


C"0(2- 1) 


21.5 ±0.02 


2.36 ± 0.06 


0.85 ± 0.06 


2.58 ±0.05 [18.4, 25.0] 


0.08 ± 0.01 


18.5 ±0.4 




H'-'C0+(3 - 2) 


21.7 ±0.08 


2.92 ± 0.53 


0.37 ± 0.05 


1.25 ±0.06 [18.5, 24.2] 


0.99 ± 0.45 


4.6 ±0.3 




DCO+(3-2) 


22.0 ±0.2 


0.56 ± 3.49 


0.04 ± 0.01 


0.04 ±0.02 [21.1, 22.8] 


0.05 ± 0.01 


4.6 ±0.3 




N2H+(3 - 2) 


21.9 ±0.02 


3.56 ± 0.02 


2.47 ± 0.27 


17.62 ±0.12 [15.6, 28.6] 


16.92 ±7.28 


8.2 ±0.1 




N2D+(3 - 2) 


22.0 ±0.1 


0.80 ± 0.20 


0.07 ±0.01 


0.07 ±0.01 [21.5, 23.0] 


0.06 


8.2 ±0.1'' 




C-C3H2(33.0~22.i) 


21.3±0.1 


1.10 ±0.26 


0.08 ±0.01 


0.07 ± 0.02 [20.6, 22.4] 


0.45 ±0.19 


3.4 ±0.1 




CH3NH2(72,5 - 7,,5)' 


21.1 ±0.2 


5.07 ± 0.58 


0.05 ± 0.01 


0.24 ±0.03 [16.7, 25.1] 


0.003 


50.4 




0CS(19 - ISf 


21.7 ±0.2 


-II- 


-II- 


-II - [17.4, 25.4] 


0.01 ± 0.0003 


13.2 ±0.4 




"CS(5 -4) 


22.3 ± 0.2 


2.69 ± 0.38 


0.05 ± 0.01 


0.18 ± 0.02 [19.7, 24.7] 


0.07 ± 0.01 


4.7 ± 0.2 




SiO(6 - 5) 


21.5 ±0.4 


6.66 ± 1.28 


0. 1 1 ± 0.04 


1.30 ±0.11 [16.5, 25.3] 


0.42 ± 0.27 


4.0 ±0.7 


G015.05+00.07 MMl 


C"0(2 - 1) 


24.6 ± 0.03 


2.22 ± 0.24 


0.74 ± 0.04 


2.05 ± 0.05 [21.6,26.9] 


0.08 ± 0.01 


16.4 ± 1.7 




H''C0+(3 - 2) 


24.5 ± 0.1 


1.66 ± 0.24 


0.08 ± 0.01 


0.20 ± 0.02 [22.8, 26.8] 


0.11 ± 0.01 


4.9 ± 0.1 




DC0+(3 - 2) 


24.1 ± 0.1 


0.67 ± 0.15 


0.08 ± 0.003 


0.07 ±0.01 [23.3, 25.2] 


0.09 ± 0.01 


4.9 ± 0.1 




N2H+(3 - 2) 


24.7 ± 0.1 


4.58 ± 0.29 


0.48 ± 0.06 


2.47 ± 0.09 [19.9, 29.3] 


0.98 ± 0.14 


5.2 ± 0.1 




N2D+(3 - 2) 


24.1 ± 0.1 


0.53 ± 0.13 


0.06 ± 0.01 


0.03 ± 0.01 [23.4, 25.4] 


0.12 


5.2 ±0.1'' 


IRDC 18151-1208 MM2 


C"0(2 - 1) 


29.5 ± 0.02 


2.75 ± 0.15 


1.03 ± 0.04 


3.33 ± 0.05 [26.5, 32.8] 


0.09 ± 0.01 


20.8 ± 1.9 




H''C0*(3 - 2) 


29.5 ± 0.04 


2.47 ±0.11 


0.64 ± 0.05 


1.63 ± 0.06 [27.4, 31.8] 


0.44 ± 0.10 


6.8 ± 0.4 




DCO+(3-2) 


29.4 ± 0.03 


1.43 ±0.10 


0.33 ± 0.03 


0.75 ±0.03 [27.3, 31.3] 


0.15 ± 0.03 


7.2 ± 0.6 




C-C3H2(33,o-22.i) 


29.5 ± 0.2 


2.87 ± 0.49 


0.07 ± 0.02 


0.18 ± 0.05 [27.2, 31.5] 


0.08 ± 0.02 


4.8 ± 0.3 




N2D+(3 - 2) 


29.4 ± 0.1 


1.03±0.18 


0.11 ± 0.03 


0.27 ± 0.03 [28.3, 30.9] 


0.04 


14.0 ± 2.0' 




C2D(3 - 2)« 


29.6 ±0.3 


2.68 ± 0.77 


0.07 ± 0.02 


0.45 ± 0.07* 


0.02 


20.8 




CH2CHC"N(252.24 - 250.25)* 


30.1 ±0.5 


3.15 ±0.97 


0.04 ± 0.02 


0.14 ±0.01 [27.8, 33.3] 


0.001 


101.2 




U216452.511 


29.4 ±0.2 


2.06 ± 0.49 


0.05 ± 0.02 


0.13 ± 0.02 [28.0, 31.3] 








SiO(6 - 5) 


30.5 ± I.O 


46.80 ±2.36 


0.09 ± 0.03 


1.93 ±0.20 [-6.1, 72.6] 


0.09 ± 0.01 


5.7 ± 0.1 


G015. 31-00.16 MM3 


C"0(2 - 1) 


31.0 ±0.03 


0.80 ± 0.09 


0.41 ± 0.07 


0.80 ± 0.04 [28.9, 32.3] 


0.10 ± 0.01 


10.2 ± 0.6 




H''CO+(3 - 2) 


31.4 ±0.2 


0.89 ±0.14 


0.03 ± 0.01 


0.05 ± 0.02 [30.2, 32.5] 


0.18 ± 0.01 


3.5 ± 0.1 




DC0+(3 - 2) 


30.8 ±0.2 


0.87 ± 0.15 


0.03 ± 0.01 


0.06 ± 0.01 [29.8, 31.9] 


0.13 ± 0.02 


3.4 ±0.1 




N2H+(3 - 2) 


30.9 ± 0.04 


0.86 ± 0.18 


0.26 ± 0.02 


0.31 ± 0.02 [29.9, 31.9]' 


4.64 ± 1.30 


3.9 ± 0.1 


IRDC 18182-1433 MM2 


C"0(2- 1) 


41.0 ±0.02 


1.32 ±0.08 


0.74 ± 0.08 


1.82 ±0.03 [38.0, 43.0] 


0.13 ±0.01 


12.6 ±0.5 




H''CO+(3 - 2) 


40.7 ± 0.1 


1.20 ± 0.34 


0.11 ± 0.02 


0.22 ± 0.02 [39.4, 42.5] 


0.50 ± 0.03 


3.8 ±0.1 




DCO+(3-2) 


41.1 ± 0.3 


2.08 ± 0.90 


0.09 ± 0.01 


0.20 ±0.03 [38.6, 43.3] 


0.28 ± 0.02 


3.8 ± 0.1 




C-C3H2(33,0 - 22,l) 


40.9 ± 0.1 


1.58 ± 0.31 


0.07 ± 0.01 


0.12 ± 0.02 [39.3, 42.4] 


1.29 ± 0.13 


3.1 ± 0.1 




CH3CUCH3(09.7 - 148.6) ■ 


41.3 ±0.1 


1 OA 1 A 1 /I 

1.80 ± 0.14 


A 1 -7 1 A Al 
0. 1 / ± 0.01 


A "3 c 1 A Al rin c /ii ni 
0.35 ± 0.03 [39.5, 42.9] 


0.01 


63.9 


IRDC 18223-1243 MM3 


C"0(2 - 1) 


45.5 ±0.04 


3.09 ±0.10 


0.56 ± 0.04 


1.97 ±0.05 [42.1, 49.0] 


0.06 ± 0.01 


18.2± 1.1 




H"CO+(3 - 2) 


45.5 ±0.04 


1.81 ±0.26 


0.39 ± 0.03 


0.83 ± 0.04 [43.4, 47.0] 


0.47 ± 0.05 


5.5 ±0.1 




DCO+(3-2) 


45.9 ±0.1 


1.83 ±0.18 


0.11 ±0.01 


0.23 ± 0.02 [43.4, 47.6] 


0.08 ± 0.01 


5.6 ±0.2 




N,D+(3-2) 


45.3 ±0.3 


1.13 ±0.63 


0.08 ± 0.02 


0.20 ± 0.03 [44.2, 47.7] 


0.05 


9.4 ± 0.3,^ 




C-C3H2(33,o-22,,) 


45.7 ±0.1 


1.84 ±0.15 


0.15 ±0.03 


0.22 ± 0.02 [43.6, 47.0] 


0.38 ± 0.05 


4.0 ±0.1 




SiO(6 - 5) 


44.4 ± 0.6 


8.22 ± 1.26 


0.07 ± 0.03 


0.84 ±0.08 [38.0, 52.6] 


0.11 ±0.01 


4.9 ±0.1 


ISOSS J 18364-0221 SMMl 


C"0(2- 1) 


34.8 ± 0.02 


1.60 ±0.05 


0.79 ± 0.03 


1.69 ±0.03 [32.5, 37.0] 


0.17 ±0.01 


10.9 ± 0.2 




H'3CO+(3 - 2) 


34.7 ± 0.05 


0.86 ±0.14 


0.30 ± 0.04 


0.44 ±0.04 [33.5, 36.1] 


1.07 ±0.38 


4.3 ± 0.2 




DC0+(3 - 2) 


34.9 ± 0.03 


1.14±0.10 


0.33 ± 0.02 


0.50 ±0.02 [33.5, 36.3] 


0.69 ± 0.20 


4.4 ±0.2 



Notes. Intensities are integrated over the velocity range given in square-brackets. In the cases of 118102 MMl, 118151 MM2, and 118223 
MM3, 92.5% of the total N2D+(3 - 2) hf component's intensity lie within the quoted velocity range. For G015.05 MMl the fraction is 81.7%. 
**' For H'^'CO*, DCO*, NaH^, and N^D^ To is the optical thickness in the centre of a hypothetical unsplit line. For C2D tq is the sum 

of the peak optical thicknesses of all the hf components. See Sect. 4.2 for details on Tex determination. T^,, is assumed to be the same 
as for N2H*(3 - 2). Another candidate for this spectral line is the 0CS(19 - 18) transition. Alternatively, the two lines could be blended. 
It could also be due to C"*0(2 - 1) seen in the image band. TexLNzH+Cl - 0)] from SSK08. Another detected hf group of C2D(3 - 2) 
is blended with CH2CHC'^N(252,24 - 25o,25). Integrated intensity over the two detected hf groups between [27.7, 30.7] and [-49.5, -46.5] 
km s"'. *'* 92.6% of the hf component's intensity lie within the quoted velocity range. This line may also be due to C'*0(2 - 1) seen in the 
image band. 

molecule (2.8 for the He/H abundance ratio of 0.1), and nin is 
the mass of a hydrogen atom. 

The values were determined from the (sub)mm maps 
shown as contours in Fig.[T] The angular resolutions of the maps 
are different: the beam size of the SCUBA 850 fj.m, Bolocam 1 . 1 
mm, and MAMBO 1.2 mm data are 14'.'4, 31", and U", respec- 
tively. The values of 5 y were taken from the literature as follows: 
forI18102 MMl, 118151 MM2, 118182 MM2, and 118223 MM3 
we used the MAMBO 1.2-mm flux densities from Beuther et 
al. (2002a); for G015.05 MMl and G015.31 MM3 we used the 
MAMBO-II 1 .2-mm flux densities from Rathborne et al. (2006); 
for J18364 SMMl we used the SCUBA 850-//m flux density 
from Birkmann et al. (2006). These flux densities are given in 
Col. (4) of Table H 

As the dust temperature of the clumps, we used the gas ki- 
netic temperatures listed in Col. (2) of Table|5] and assumed that 
T^dust - T^kin- We extrapolated the values of Ky from the Ossenkopf 
& Henning (1994) model for dust grains with thin ice mantles. 



4.1.3. Clump masses, radii, and H2 column and number 
densities 

We calculated the masses and beam-averaged H2 column densi- 
ties of the clumps using the formulas 

M=^^, (1) 

y^dust 

mUi) = . (2) 

In the above formulas, /jj"'^* and S y are the peak surface bright- 
ness and integrated flux density of the (sub)mm dust emission; 
note that /jj"^' = ^^'^''''/Qbeam, where Sl'''^^ is the peak flux den- 
sity and Qbeam is the solid angle of the telescope beam. Bv(rdu.st) 
is the Planck function for a dust temperature Tdust, Ky is the dust 
opacity per unit dust mass, Rd = (Mdust /A/gas) is the average 
dust-to-gas mass ratio, //h, is the mean molecular weight per H2 
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coagulated for 10 yr at a gas density of hh - n(H) + 2n(H2) - 
2n(H2) = 10*" cm"^ [their Table 1, Col. (6)]. This is expected to 
be a reasonable dust model for the sources within cold and dense 
IRDCs. At the wavelengths considered in the present work, these 
values are /cgso/zm - 0.197 m^ kg"', /ci.imm = 0.121 m^ kg"', 
and Ki2mm = 0.106 m^ kg"'. The canonical value 1/100 was 
adopted for R^- We note that the dust opacities are likely to 



be uncertain by a factor of > 2 ( ,Ossenkopf & Henning 1994 
IMotte & Andre 20011 1. For comparison, Williams et al. (2004), 
Enoch et al. (2006), and Rygl et al. (2010) used the same dust 
model as here and interpolated the values of Ky at 850 fim, 1 . 1 
mm, and 1.2 mm to be 0.154, 0.114, and 0.1 m^ kg"', respec- 
tively. These are slightly smaller than our values, and the dif- 
ference is likely to be caused by different extrapolation methods 
(e.g., log-interpolation) and/or different dust emissivity index, 
used by the authors to determine the value of a:v oc v'^. 

In the papers by Beuther et al. (2002a) and Rathborne 
et al. (2006), the reported clump sizes refer to their FWHM 
sizes (diameters) resulting from two-dimensional Gaussian fits. 
However, the integrated flux densities used to calculate the 
masses refer to larger clump areas, A. Thus, in order to properly 
calculate the volume-average H2 number density, (n(H2)), one 
way is to use the so-called efl'ective radius of the clump, /?eff - 
yfAjji (see, e.g., Miettinen & Harju 2010) 1. An alternative way to 
calculate (n(H2)) is to compute the amount of mass within the 
FWHM contour and use the FWHM radius. For sources with a 
Gaussian shape, the mass within the FWHM contour is a fraction 
In 2 ^ 0.693 of the total mass (see I Kauffmann & Pillai 20T0b . 
Following Kauffmann & Pillai (2010), we reduce the clump 
masses by the above factor and use half the FWHM size as the 
effective radius of the clump. The density (n(H2)) was then cal- 
culated using the formula 



<«(H2)> 



ip) 



(3) 



where (p) = M/ (a/SttR^^^^ is the mass density, and M refers 
to the reduced clump mass as described above. For the source 
J18364 SMMl, we scaled the effective radius » 0.2 pc reported 
by Birkmann et al. (2006) using the revised distance (/?eif = 0.23 
pc), and used the total mass within the clump area. 

The results of the above calculations are presented in 
Cols. (5)-(8) of Table|5] The uncertainties in the derived quanti- 
ties were propagated from the uncertainties in d and Tkin, but we 
note that the uncertainty in (sub)mm dust opacity is the largest 
source of error in M and A^(H2). 

The derived values of M, N(H2), and (n(H2)) are mostly 
within a factor of < 2 of those derived by Beuther et al. (2002a; 
where the revised equations by Beuther et al. (2005b) are em- 
ployed), Rathborne et al. (2006), Birkmann et al. (2006), and 
Hennemann et al. (2009). We note that besides the different 
source distances, the dust parameters used by the authors in the 
above reference studies were different than here. Beuther et al. 
(2002a, 2005b) used Tdust values (35-50 K) resulting from the 
cold part of the source SEDs; the SEDs were constructed from 
the IRAS and mm data of the main core in the source region 
which explains the rather high dust temperatures. In addition, 
Beuther et al. (2002a, 2005b) used the grain radius, mass den- 
sity, and emissivity index values of a = 0. 1 yum, pdust = 3 g cm"-', 
and 13-2, respectively, in the calculation of mass and column 
density. Following Hildebrand (1983), these values correspond 
to K\2mm ~ 0.5 m^ kg"'. Rathborne et al. (2006) used the values 
T^dust = 15 K and K\2mm - 0.1 m^ kg"'. For the source J18364 
SMMl, Birkmann et al. (2006) used Tdust = 16.5+3„ K instead 



of Tkin, and they had /csso/im = 0.180 m kg ' which is slightly 
smaller than ours even the adopted dust model was the same. 

4.2. Molecular column densities and fractional abundances 

The beam-averaged column densities of C'^O, H'^CO^, DCO^, 
N2H+, SiO, OCS, '""CS, and 0-C-C3H2 were derived using a one- 
dimensional spherically symmetric non-LTE radiative transfer 
code called RADEX (see Sect. 3.2). RADEX uses the method of 
mean escape probability for an isothermal and homogeneous 
medium. The molecular data files (coUisional rates) used in 
the RADEX excitation analysis were taken from the LAMDA 
database (ISchoier et al. 20051) . The C'^O, H'^^CO+, DCO+, and 
N2H^ transitions are treated as a hypothetical unsplit transition. 
The input parameters in the off-line mode of RADEX are the gas 
kinetic temperature, H2 number density, and the width (FWHM) 
and intensity of the spectral line. We used the values of Tkin and 
(n(H2)) listed in Table |5] However, we multiplied the densities 
by 1 .2 (He/H2 = 0.2) to take the collisions with He into account 
(see Sect. 4. 1 of the RADEX manuaQ). As the input line intensity 
we used the main-beam brightness temperature Tmb - T^/rjMB- 
When the source is resolved, Tmb is equal to the Rayleigh- Jeans 
equivalent radiation temperature, Tr. The simulations aim to re- 
produce the observed line intensity and yield the values of the 
line peak optical thickness (tq), excitation temperature (Tex), and 
the total column density of the molecule (A^tot)- We varied T]^i„ 
and {n(H2)) according to their errors to estimate the uncertain- 
ties associated with tq, T^^, and A^tot- The lines appear to be op- 
tically thin in most cases. The optical thickness is tq > 1 for all 
the N2H+ fines, for H'3C0+ towai'ds 118102 MMl and J18364 
SMMl, and for C-C3H2 of 118182 MM2. Tfie Tex values for 
C'^O are close to T^„, indicating tfiat tfie lines are nearly tfier- 
malised. Also, rex(H''CO+) is found to be equal to rex(DCO+) 
witfiin tfie errors. We note tfiat tfie column density determination 
is in some cases fiampered by self-absorbed line profiles. 

For tfie rest of tfie observed molecules, tfiere are no molecu- 
lar data files available in tfie LAMDA database. Tfie line op- 
tical tfiicknesses and total beam-averaged column densities of 
tfiese molecules were determined tfirougfi LTE modelling witfi 
CLASS/Weeds. Tfie input parameters for a Weeds model are 
Ntot, T^ex, source size (©s), linewidtfi (FWHM), and offset from 
tfie reference-cfiannel velocity. Tfie linewidtfi is directly deter- 
mined from tfie observed line profile, so tfiere are basically tfiree 
free parameters left (Mot, ^ex, 6s)- Some of tfie model parame- 
ters may be degenerate, and cannot be determined independently 
(IScfiilke et al. 20061 IMaret et al. 201 11 ). Tfie source size is de- 
generate witfi excitation temperature in tfie case of completely 
optically tfiick lines (t » 1), and witfi column density if tfie 
lines are completely optically tfiin (r 1). We assumed tfiat 
tfie source fills tfie telescope beam, i.e., tfiat tfie source solid 
angle, D.^ oc 6^, equals tfie beam solid angle, Q.^ ^hpbw '^^ 
da - ^^HPBW- Wfien tfie beam filling factor is unity, tfie line brigfit- 
ness temperature is - Tmb [see Eqs. (1) and (2) in Maret et al. 
(201 1)]. For N2D^, we used as Tex tfie values obtained for N2H^ 
from RADEX simulations or from SSK08 (Sect. 3.2). For tfie rest 
of tfie transitions we adopted tfie value Tex - E^/ks for linear 
rotors, and Tex = 2/3 x E^lk-^ for asymmetric top rotors. Tfie 
adopted rotational excitation temperatures give a lower limit to 
tfie total beam-averaged column densities CHatcfiell et al. 1998] l. 
Tfie input A^tot was tfien varied until a reasonable fit to tfie line 
was obtained (see Fig. |7]l. We note tfiat for tfie complex or- 
ganics we fiave detected only one transition per species. We 
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Table 5. Physical properties of the sources. 



Source 


^kin' 




7 b 

J dust 


s/ 


M 








<n(H2)> 




[K] 




[K] 


[Jy] 


[Mq] 


[10^2 


cm ^] 


[pc] 


[10^ 


cm ^] 


IRDC 18102-1800 MMl 


21.3 ± 


1.6 




3.3 


357 ± 


163 


6.5 


tO.7 


0.34 


2.9 


± 1.3 


G015.05+00.07 MMl 


17.2 ± 


2.1 


11.0-36.0 


0.47 


63 + 


31 


I.O 


±0.2 


0.15 


5.9 


±2.0 


IRDC 18151-1208 MM2 


21.1 ± 


2.0 




2.6 


285 ± 


111 


12.9 


± 1.6 


0.17 


18.5 


±7.2 


G015.31-00.16 MM3 


13.7 ± 


2.8 




0.74 


183 + 


83 


0.8 


tO.3 


0.38 


1.1 


±0.5 


IRDC 18182-1433 MM2 


15.3 ± 


1.5 




0.3 


86 ± 


23 


5.2 


±0.8 


0.26 


1.6 


±0.4 


IRDC 18223-1243 MM3 


18.7 + 


1.3 


15.0/18.0 


0.8 


173 + 


43 


1.4 


tO.I 


0.18 


9.5 


±2.3 


ISOSS J18364-0221 SMMl 


11.4 


16 5+'' '' 


2.11 


169 + 


54 


7 


.6 


0.23 


6.4 


±2.0 



Notes. Calculated from r„,(NH3) from Krause (2003) and SSK08. Krause (2003) reported a slightly higher Ti^„ value of 1 1.8 K for J18364 
SMMl. <*' For G015.05 MMl, is from Rathborne et al. (2010). For 118223 MM3, the two Tj^^ values, 15 and 18 K, are from Beuther & 
Steinacker (2007) and Beuther et al. (2010), respectively. The T^^^t value for J18364 SMMl is from Birkmann et al. (2006). For those sources 
whose namestart with "IRDC", 5, refers to MAMBO 1 .2-mm flux d ensity from B euther et al. (2002a). ForG015.05 MMl andG015.31 MM3, 
the quoted value is the MAMBO-II I.2-mm flux density ( Rathborne et al. 2006t . In the case of J18364 SMMl, we give the SCUBA 850-yum 
flux density from Birkmann et al. (2006). 



cannot thus apply the method of rotation diagram to derive 
rotational temperatures and molecular column densities (e.g.. 
Goldsmith & Langer 1999] ). 

We also determined the HCO^ column density from the co- 
lumn density of H'^CO^. For this calculation, it was assumed 
that the carbon-isotope ratio ['^C]/['-'C] depends on Rqc accor- 
ding to the relationship given by Wilson & Rood (1994): 



=7.5x/?Gc[kpc]+7.6. (4) 

For the Rcc values considered here, 5.1 -6.3 kpc, the above ratio 
lies in the range ~ 46 - 550 

For 118151 MM2, 118182 MM2, and 118223 MM3 we do 
not have N2H^ data. For these sources we computed the N2H^ 
column density from the 7=1-0 line parameters (Tex and 
t) from SSK08 [see, e.g., Eq. (10) in Miettinen et al. (2010)]. 
We obtain about 1.1-1.2 times higher N2H^ column densities 
compared to SSK08 who used the optically thin approximation 
(see Sect. 5.1 for a further discussion). 

We calculated the fractional abundances of the molecules by 
dividing the molecular column density by the H2 column den- 
sity: x(mol) = A^(mol)/A^(H2). For this purpose, the values of 
A^(H2) were derived from the (sub)mm dust continuum maps 
smoothed to the corresponding resolution of the line observa- 
tions. With a resolution of 31", the Bolocam 1.1 -mm data could 
not be smoothed to correspond the resolution of the line obser- 
vations. In these cases we used the original Bolocam data; the 
31" resolution is in most cases comparable to that of the line ob- 
servations (22"3 - 28"9), and thus we do not expect this to be 
a significant source of eiTor. The derived column densities and 
abundances are listed in Tables |6] and [T] The abundance errors 
were derived by propagating the errors in A^(mol) and A^(H2). 



IRDC •8182-1433^^^ 









2,31085 10^ 2,3108 10' 2,3-075 10^ 

Rest Frequency (MHz) 



Fig. 7. Example of the Weeds LTE modelling outlined in 
Sect. 4.2. The synthetic model spectrum is overlaid as a red line. 
Note that the intensity scale is Tmb for modelling purposes. 



Table 7. Column densities and fractional abundances of the other 
candidate species observed. 



Source 


Molecule 


Nm [cm--] 


.V = A',o,/A'(H2) 


IRDC 18102-1800 MMl 


0-C-C3H2 


1.7 ± 0.8 X 10'-' 


1.8 ± 0.8 X lO-' 




CH3NH2° 


5.0 X 10'*'' 


5.1 ±0.6 X 10-* 




DCS" 


2.3±0.8x 10'5 


2.4 ± 0.9 X 10-' 




'^CS 


6.0 ±2.9x10'' 


6.2 ±3.1x10-' 




SiO 


1.2 ±0.8x10" 


1.1 ±0.8 X 10-' 


IRDC 18151-1208 MM2 


0-C-C3H2 


7.0 ± 1.7 X 10'- 


1.9 ±0.5 X 10-'" 




C2D' 


1.5 X 10'" 


4.2 ±0.5 X 10-'" 




CH2CHC''^N 


2.5 X lO'"* 


6.9 ± 0.9 X 10-' 




SiO 


1.1 ±0.3x 10'^ 


2.5 ± 0.8 X 10-* 


IRDC 18182-1433 MM2 


0-C-C3H2 


7.7 ± 0.8 X 10" 


1.2 ± 0.2 X 10-* 




CH3COCH3 


2.0 X 10'^'' 


4.1 ±0.6 X 10-' 


IRDC 18223-1243 MM3 


0-C-C3H2 


2.2 ± 0.3 X 10'-' 


1.6 ± 0.3 X 10-' 




SiO 


4.0±0.6x lO'-t 


2.9 ± 0.5 X 10-* 



Notes. The corresponding line is blended with OCS. * ' Calculated 
by assuming that Tex = lI'ixE^jk^.''^'' The other C2D hf group detected 
is possibly blended with CH2CHC'^N. Calculated by assuming that 



Spectral lines of ''C-isotopologue of HCO^ are likely to be opti- 
cally thick. Therefore, the HCO^ deuteration can be better investigated 
through the DCO^/H'^^CO^ column density ratio. However, a caveat 
should be noted here. The HCO^ molecules are produced directly from 
CO (see reactions 3 and 7 in Table [Sj. On the other hand, at low tem- 
perature, CO is susceptible to the exothermic isotopic charge exchange 
rea ction "C+ + '^CO ^ ^^C^ + '^CO + AE, where AE/ks = 35 
K dWatson et al. 1976] >. This is expected to cause considerable '^C- 
fractionation in cold and dense gas, which complicates the deuteration 
analysis. 



4.3. CO depletion and deuterium fractionation 

To estimate the amount of CO depletion in the clumps, we cal- 
culated the CO depletion factor, /d, following the analysis pre- 
sented in the paper by Fontani et al. (2006). If x(CO)can is the 
"canonical" (undepleted) abundance, and x(CO)obs is the ob- 
served CO abundance, fu is given by 

_ x(CO)can 

JD - ■ (->) 

x(CO)obs 
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Table 6. Molecular column densities and fractional abundances with respect to H2. 



Source 


A'(C"0) 


A'(N2H+) 


A'(N2D+) 


A'(H"CO+) 


iV(HCO+) 


A'(DCO+) 




[10'^ cm--] 


[10'^ cm-^] 


[10" cm-2] 


[10'^ cm-2] 


[10'* cm--] 


[10" cm--] 


IRDC 18102-1800 MMl 


1.4 ± 


0.1 


30.1 ± 13.3 


4.9 


1.9 ± 1.0 


9.7 ±5.1 


1.9 ±0.9 


G015.05+00.07 MMl 


1.2 ± 


0.1 


3.3 ±0.6 


9.2 


0.1 ±0.03 


0.5 ±0.2 


3.0 ±0.6 


IRDC 18151-1208 MM2 


2.1 ± 


0.1 


5.0 ± 1.1" 


5.2 


0.5 ±0.1 


2.6 ±0.5 


8.2 ± 1.9 


G015.31-00.16 MM3 


0.3 ± 


0.1 


3.8 ±0.6 




0.3 ± 0.04 


1.5 ±0.2 


11.0±2.1 


IRDC 18182-1433 MM2 


0.7 ± 


0.1 


2.1 ±0.4" 




0.7 ±0.1 


3.2 ±0.5 


45.4 ±5.4 


IRDC 18223-1243 MM3 


1.2 + 


0.1 


5.1 ±0.5" 


6.7 


0.4 ±0.1 


1.9 ±0.5 


6.6 ± 0.9 


ISOSS J18364-0221 SMMl 


1.1 ± 


0.1 






0.6 ± 0.2 


3.3 ± 1.1 


39.2 ± 13.8 




x(C"0) 


x(N2H+) 


x(N2D+) 


x(H'^CO+) 


x(HCO+) 


x(DCO+) 




[10- 


'] 


[10-'] 


[10-"] 


[io-'«] 


[10-'*] 


[10-"] 


IRDC 18102-1800 MMl 


1.5 ± 


0.2 


27.7 ± 12.6 


5.0 ±0.5 


18.2 ±9.8 


9.3 ±5.0 


2.0 ± 1.0 


G015.05-i-00.07 MMl 


1.2 ± 


0.2 


3.4 ±0.9 


9.5 ± 1.7 


1.0 ±0.4 


0.5 ±0.2 


3.1 ±0.8 


IRDC 18151-1208 MM2 


0.6 ± 


0.1 


0.9 ±0.2 


1.3 ±0.2 


1.2 ±0.3 


0.6 ±0.1 


2.3 ±0.6 


G015.31-00.16 MM3 


0.4 ± 


0.2 


4.5 ± 1.6 




3.5 ± 1.2 


1.8 ±0.6 


13.0 ±4.8 


IRDC 18182-1433 MM2 


1.1 ± 


0.2 


2.5 ± 0.6 




9.9 ±2.1 


4.5 ± 1.0 


72.8 ± 14.6 


IRDC 18223-1243 MM3 


0.9 ± 


0.1 


3.7 ±0.5* 


4.9 ±0.5 


2.9 ±0.8 


1.4 ±0.4 


4.8 ±0.8 


ISOSS J18364-0221 SMMl 


1.6 ± 


0.1 






7.3 ±2.4 


4.0 ± 1.3 


58.7 ± 20.7 



Notes. Calculated from the SSK08 parameters (18" resolution). The corresponding H2 column density was estimated from the Bolocam 
map of 31" resolution. 



The "canonical" CO abundance at the galactocentric dis- 
tance Rqc was calculated using the relationship [Eq. (7) in 
IFontanietar2006J 



jc(CO)can = 9.5 X IQ-^e' 



(6) 



Note that this relationship results from using the value Rq = 8.5 
kpc, whereas our Rqc values are computed using Rq - 8.4 kpc. 
Such a small discrepancy is however negligible. At Rcc - 8.5 
kpc, the above relationship gives the standard value 9.5x1 0^^ for 
the abundance of the main CO isotopologue in the solar neigh- 
bourhood ( [Frerking et al. 1982| . To calculate the "canonical" 
C'^O abundance we take into account that the oxygen-isotopic 
ratio [^^0]/['^0] depends on Rqc according to the relationship 
dWilson & RoodT994j 



ri6, 



O] 



[ISO] 



= 58.8 x/?Gc[kpc]H- 37.1. 



(7) 



For the Rqc values of our clumps, the above ratio range from 
about 337 to 407.5. When the above relationship is combined 
with the ['^0]/['^0] ratio, for which we use the standard value 
3.52 ([Frerking et al. 1982[), the value of x(C'^0)can can be cal- 



culated as 



ji:(C"0)c 



x(CO)c 



[i80]/[i70] X [i60]/[i80] 

x(CO)can 

3.52 X (58.8 x/?GC [kpc] H- 37.1) 



(8) 



The depletion factor fu is then calculated from fu - 
4C'^0)can/x(Ci''0)obs, and the results, /d = 0.6±0.1 -2.7± 1.3, 
are listed in Col. (2) of Table |9] The ±-error quoted was calcu- 
lated by propagating the uncertainty in x(C'^0)ob,s- The low va- 
lues of /d indicate that CO is not significantly depleted, if at all, 
in our clumps. 

The degree of deuterium fractionation in HCO^ and 
N2H^ was calculated by dividing the column density of the 
deuterated isotopologue by its normal hydrogen-bearing form 
as /?d(HCO+) = A?(DCO+)/A^(HCO+) and /?d(N2H+) = 
A^(N2D+)/A^(N2H+). The values of /?d(HCO+) and /?d(N2H+) 



are in the range 0.0002 ± 0.0001 - 0.014 ± 0.003 and 0.002 ± 
0.001 - 0.028 + 0.005, respectively; the error in Ru was de- 
rived from the errors in the corresponding column densities [see 
Cols. (3) and (4) of Table|9l. 

4.4. lonisation degree 

We have derived fractional abundances of different ionic 
species and their different isotopologues. In principle, a lower 
limit to the ionisation degree can be obtained by sim- 
ply summing up these abundances (e.g., ICaselli et al. 20021 
llVTiettinen et al. 2009l l: 



x(e) > x(HCO+)-Hx(H'^CO+)H-x(DCO+)H-x(N2H+)-Hx(N9D+). 

(9) 

This is based on the gas quasi-neutrality: the electron abundance 
equals the difference between the total abundances of the cations 
and the anions. The resulting values, x(e) > 2; x(ions),, are 
listed in Col. (5) of Tabled These lower limits are in the range 
7xl0-''-1.2xl0-^ 

The fractional ionisation in dense molecular clouds can 
also be inferred by employing ion-molecule chemical reaction 
schemes with the observed molecular abundances. All the reac- 
tions considered here, and notes on the associated rate coeffi- 
cients, are given in Table |8] 

At first, we derive another estimate for a lower limit to x(e) 
through the ionisation balance determined by , HCO^, N2H^, 
and electrons following Qi et al. (2003). The and N2H^ are 
mainly destroyed by CO; at steady state, reactions 3, 1 1, and 15 
in Table |8] lead to the following equation for the lower hmit to 
x(e) [cf. Eq. (6) in Qi et al. (2003)]: 



x(e) > 



^i5x(N2H+)x(CO) 

y8llX(HCO+) 



(10) 



The rate coefficients k\5 and jSu were adopted from the UIVIIST 
databasJ'll Woodal l et al. 2007l l. The derived values are listed in 
Col. (6) of Tabled For 118151 MM2 and 118223 MM2 the x(e) 
values derived from Eq. (fTOl) are similar to those computed from 
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Eq. ©.For 118102 MM1,G015.31MM3, and 118182 MM2 the 
summed abundance of different ionic species is clearly higher 
(by factors ~ 3.5 - 16.3) than the lower limit to .^^(e) resulting 
from Eq. ([TOll. In the case of G015.05 MMl, on the other hand, 
the summed ionic abundance is about six times lower. These dis- 
crepancies are not surprising because the chemical scheme be- 
hind Eq. ( fTOb is certainly oversimplified. 

Next, we determine the degree of ionisation by utilising the 
abundance ratios RoiUCO^) = [DCO+]/[HCO+] and Ra = 
[HCO^]/[CO]. The first studies of fractional ionisation based on 
the above ratios were carried out over three decades ago (e.g., 
IGuelin et al. T9771 IWatson et al.T 978; Wootte n et al. 197% . A 
similar analysis was subsequently applied in the papers by 
CaselU et al. (1998), WilHams et al. (1998), Bergin et al. (1999), 
Anderson et al. (1999), and Caselli (2002). 

We note that the following analysis includes several assump- 
tions: /) HCO^ is mainly produced in the reaction between 
Hj and CO (reaction 3 in Table [8]i; ii) all deuteration is due 
to the reaction between H3 and HD (reaction 2); Hi) the pre- 
sence of atomic deuterium, which could (slightly) increase the 
deuteration degree, is ignored; iv) ionic species are destroyed 
mainly by electrons, the most important neutrals (CO and O), 
and negatively charged dust grains; v) except CO and O, we 
neglect the contribution of some other neutral species, such as 
N2, O2, and H2O, in the destruction of and HaD^; and v;) 
we ignore the effects of refractory metals (Anderson et al. 19991 
ICaselli et al. 2002l l. Concerning part v), the abundances of neu- 
trals, such as N2 and H2O, are poorly known and/or low. For 
instance, the first results from Herschel have shown that H2O 
abundance is relatively low in high-mass star-forming regions 
(~ 10'"' - 10'^ van der Tak et al. 201 0; Mar seille et al. 2010al 
IChavarrfa et al. 20101) . Bergin et al. (1999) varied the nitrogen 
abundance in their chemical model and found that it does not af- 
fect the electron abundance. Also, the destructive reaction with 
O2 is very slow {k = 9.3 x lO""' cm^ s"'; UMIST). 

By writing steady-state equations for the abundances of 
H2D+ (reactions 2 and 7-10), DCO+ (reactions 7 and 13-14), 
and HCO^ (reactions 3, 7, and 1 1-12), it can be shown that 



7Jd(HCO+) = \, 



fc+x(HD) 



3 kjxiCO) + p^xio) + kgxiO) + kwx(g) 
Solving x(e) from the above formula yields 



(11) 



x(e) 



1 



k+x(liD) 



3/?d(HCO+) 



• hx(CO) - kgx(0) - kiox(g) 



(12) 



This represents the upper limit to ji:(e). We have assumed that 
kj - 1/3 X ki, because H2D^ can transfer a proton to CO 
producing HCO^; the rate for this is two times higher than 
for the channel producing DCO^. For the reaction rate of the 
deuteration reaction -1- HD {k+) we used the recent results 
by Hugo et al. (2009; their Table Vlll). For example, at 10 K 
the deuteration of H^ proceeds about 4.4 times faster according 
to Hugo et al. (2009) than suggested by laboratory measure- 
ments of Gerlich et al. (2002). Such a high rate was also used 
in the studies cited above, i.e., k+ = 1.5 x 10"^ cm^ s"' (e.g., 
ICaselli et al. 1998T l. We take the HD abundance to be twice the 
elemental D/H-ratio, i.e., ^HD) = 2 x [D]/[H] ~ 3 x 10"^ (e.g. 



Linsky et al. 20061 IProdanovic et al. 20101 ). Because CO is not 



found to be depleted in our clumps, we assume that this is also 
the case for atomic oxygen and use the "standard" abundance 
relative to H2 of x(0) = 3.52 x 10""^, i.e., comparable to x(CO) 



(see lCaselli et al. 19981 ICaselli 2002T l. For the grain abundance, 
jc(g), we use the value 2.64 x 10"^^, which is based on the va- 
lues fl = 0.1 yum, pgrain = 3 g cm""', and = 1/100 [see, e.g., 
Eq. (15) in Pagani et al. (2009a)]. The derived upper limits to 
x(e) are shown in Col. (7) of Tabled The values lie in the range 

~ 2 X 10"^ 3 X 10""*, i.e., much larger than the lower limits 

estimated above. 

By deriving a steady-state equation for the H3 abundance 
(see reactions 1-6 in Table[8]), and applying it in the correspon- 
ding equation for HCO^, it can be shown that 



R» = 



[^H,/«(H2)]^3 



l/^4x(e) + k3x(CO) + k^xiO) + kexig)] {j3nx(e) + knx(g)] ' 

(13) 

After Hj has formed via cosmic-ray ionisation of 
H2 (reaction 1), it quickly reacts with H2 to form H^ 
(Solomon & Werner 19711 ). Thus the H3 abundance is governed 
by the rate ^Hj- When the fractional ionisation in the source is 
determined, Eq. ( fT3T l can be used to infer the cosmic-ray ioni- 
sation rate of H2. To calculate (h^, for each source we adopted 
as x(e) the summed abundance of ionic species. Using the de- 
rived upper limits to x(e) would result in unrealistically high 
values of (u2- This suggests that the upper x(e) limits are clearly 
higher than the true values. The obtained results are shown in 
the last column of Table |9] In most cases, the values lie in 
the range ~ 1 x 10"'^ - 5 x 10"'^ s"'. However, towards 118102 
MMl we obtain a very high rate of ~ 1 x lO"'*" s"'. Because 
118102 MMl is associated with a bright MIR point-source, and 
shows no signs of CO depletion, it is probably a rapidly evolving 
source. Therefore, a steady state assumption used in the above 
analysis may be invalid for 118102 MMl. 



5. Discussion 

In the following subsections, we discuss our results, and com- 
pare them with the results from previous studies. 



5.1. Molecular column densities and abundances 

Some of the molecular column densities and fractional abun- 
dances derived here have been determined in previous studies. 
Assuming LTE conditions and optically thin emission, SSK08 
estimated N2H^ column densities from the 7=1-0 transition 
for 118102 MMl, G015.05 MMl, and G015.31 MM3 which are 
clearly lower than our values (by factors ~ 2-6). SSK08 as- 
sumed that rex(N2H+) = rrot(NH3) even though the rex(N2H+) 
values are considerably lower than r,ot(NH3) (their Table 4). On 
the other hand, rex[N2H+(l - 0)] from SSK08 ai-e very similai- 
to those from our RADEX simulations for the 7 = 3-2 transition. 
Also, the line optical thicknesses SSK08 list in their Table 4 in- 
dicate optically thick emission (t ~ 2 - 1 1). This can explain the 
discrepancy in the derived column densities. 

Mai-seille et al. (2008) modelled the N2H+, HCO+, and 
H'^CO^ emission of 118151 MM2, and obtained abundances 
which are comparable within the errors to those we have derived. 
This strengthens the reliability of our HCO^ abundances derived 
fromH'^CO^ by utiHsing the ['^C]/['^C] abundance ratio. 

Sakai et al. (2010; hereafter SSHIO) derived H1^C0+ co- 
lumn densities from the 7=1-0 transition for 118102 MMl, 
118151 MM2, and 118223 MM3. Their results ai-e otherwise sim- 
ilar to ours except that for 118102 MMl the column densities 
differ by a factor of ~ 4. SSHIO assumed optically thin emis- 
sion and that T^^{H^^CO^) = ^^((CHjOH) ± 5 K= 16.7 ± 5 
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Table 9. Depletion-, deuteration-, and ionisation parameters of the clumps. 



Source 


/d 




^d(N2H+) 


2, .v(ions), 


x(e)i 


Jc(e)„ 


in2 












[10-'] 


[10-'] 


[10-5] 


[10-" s-'] 


IRDC 18102-1800 MMl 


0.7 ± 


0.1 


0.0002 ± 0.0001 


0.002 ± 0.001 


12.3 + 6.4 


3.5 


29.2 


115 


G015.05+00.07 MMl 


0.8 ± 


0.1 


0.006 ± 0.003 


0.028 ± 0.005 


0.9 + 0.3 


5.6 


0.8 


1.2 


IRDC 18151-1208 MM2 


1.6 ± 


0.3 


0.003 ± 0.001 


0.010 ± 0.002" 


0.7 + 0.1 


0.7 


1.8 


5.0 


G015.31-00.16MM3 


2.7 + 


1.3 


0.007 ± 0.002 




2.3 + 0.8 


0.6 


0.7 


5.5 


IRDC 18182-1433 MM2 


1.1 ± 


0.2 


0.014 ± 0.003 




4.9+1.1 


0.3 


0.2 


19.7 


IRDC 18223-1243 MM3 


1.3 ± 


0.1 


0.003 ± 0.001 


0.013 + 0.001" 


1.8 + 0.5 


1.6 


1.8 


13.1 


ISOSS J18364-0221 SMMl 


0.6 ± 


0.1 


0.012 ± 0.006 




4.1 + 1.3 




0.3 


51.1 



Notes. Calculated by utilising A'(N2H+) from SSK08 (see Sect. 4.2 and 5.1). 



Table 8. Ion-molecule reactions included in the analysis of ion- 
isation degree. 



No. 


Reaction 


Note on 
rate coefRcient" 


1 


Ht + cip — ^ HJ + 
H2 + H2 — > + H 




2 


H* + HD ^ HtD+ + Hi 

^ k- 


/:+ from Hugo et al. (2009) 


3 


H+ + CO-^HCO+ +H2 


k, = lJx 10-' cm^ 


4 


H+ + e- ^ H2 + H 

H + H + H 


Pi fr om Pagani et al. (2009a) 


5 


H| +O^H20+ +H 
OH+ + H2 


ki = 1.2x 10"' cm' 


6" 


H+ + g-ig°+H2+H 

g" + H + H + H 


ku from Pagani et al. (2009a) 


7 


H2D+ +CO-^HCO+ +HD(|) 
DCO+ +H2 (i) 


k^ = l/3x/t3 


8 
9 


H2D+ +e-^H + H + D 
HD + H 
H2 +D 

H2D+ + 0H+ + HD 
OD+ + H2 


/3g fr om Pagani et al. (2009a) 
kg = ks 


10 
11 


H2D+ + g- ^ g° + H2 + D 
g° + HD + H 
g° + D + H + H 

HCO+ + e- ^ CO + H 


kio from Pagani et al. (2009a) 

^ / T-, . ,-0.69 

/;„ =2.4x10-' (3^) 


12 


HCO+ + g- ^ g° + CO + H 


ki2 from Pagani et al. (2009a) 


13 


DCO+ + e- ^ CO + D 




14 


DCO+ + g- ^ g° + CO + D 


k\i = kn 


15 


N2H+ + CO HCO+ + N2 


ki5 = 8.8 X 10-'° cm' s-' 



Notes. *"' The rate coefficients are taken from the UMIST database un- 
less otherwise stated. The temperature-dependent rates were calculated 
by using the Ti^[„ values listed in Col. (2) of Table[5]*''' The label g refers 
to the dust grains. 



K~ 3.6 X rex[Hi3CO+(3 - 2)]. SSHIO also derived SiO column 
densities from the 7 = 2-1 transition for the above three sources. 
We have obtained 36-86 times higher column densities for these 
sources from the 7 = 6-5 transition through RADEX calcula- 
tions. The spatial resolution of SiO observations by SSHIO was 
18", and as rex(SiO) they used r„t(CH30H) ± 5 K, or assumed 
that rex(SiO) = 20 K, which are over four times higher than 
rex[SiO(6 - 5)] derived here. For comparison, we performed the 
LTE modelling of SiO(6 - 5) lines in Weeds, and found that 
column densities comparable to those obtained from RADEX are 
needed to explain the observed line intensities. Of course, simi- 
larly to SSHIO, our RADEX and Weeds analyses of a single SiO 
transition include several (uncertain) assumptions. The large dis- 
crepancy in the derived A^(SiO) values could also be due to the 



IRDCs (jRagan et al. 2006 



iBeuther & Henning 2009 



fact that the 7 = 2-1 and 7 = 6-5 lines, which have completely 
different upper-state energies (6.3 K and 43.8 K, respectively), 
originate in different parts within a clump. 

In Fig. |8] we show stock charts of the derived molecular 
column densities and fractional abundances. For compa- 
rison, we plot column densities and abundances from several 
previous studies. The quoted reference studies deal with 

SSK08; IGibson et al. 20"09l 
IChen et al. 20T0l SSHIO; 
Vasyunina et al. 201 I|l, high-m ass young stellar objects 
(HlVTYSOs: ' Fontani et al. 20061 IThomas & Fuller 2008). and 
low-mass starless and protostellar cores (Butneretal. 1995}. 
The latter study was included because we found no other 
reported DCO+ surveys towards IRDCs or HMYSOs. We also 
plot values from Fontani et al. (201 1) who studied a sample of 
27 high-mass cores (including starless cores, HMYSOs, and UC 
Hii regions). The C'**0 abundances from Ragan et al. (2006) 
were converted into x(C'^O) by utilising the abundance ratio 
['^0]/[i^0]= 3.52. 

As can be seen in Fig. |8] the column densities and abun- 
dances we have derived are generally comparable to those ob- 
served in other sources. The most evident difference is the high 
SiO and '^CS column densities and abundances we have ob- 
tained compared to the IRDCs studied by SSHIO and Vasyunina 
etal. (2011). 

SSHIO found that, in general, the N(SiO)/N(H^^CO^) ratio 
is higher for the MSX dark sources than for the sources asso- 
ciated with IVIIR emission, and suggested that this ratio repre- 
sents the fraction of the shocked gas. For the sources 118102 
MMl, 118151 MM2, and 118223 MM3, SSHIO found the ratios 
A^(SiO)/A^(H'3CO+) = 2.9^\ l, 2.4+i-^, and 1.6!i-^, respectively. 
Using our column densities these ratios become very different: 
57.1 ± 46.8, 220 + 74.4, and 80 + 20, respectively. The latter 
values are in better agreement with the general trend found by 
SSHIO for a sample of 20 sources. SSHIO suggested that the 
SiO emission from the MIR dark objects originates in newly 
formed shocks, whereas the SiO emission from more evolved, 
MIR bright objects could originate in gas shocked earlier in 
time. This is supported by the fact that we see the narrowest 
SiO line towards 118102 MMl, presumably the most evolved 
source in our sample with detected SiO emission (see Sect. B.l 
of Appendix B). Follow-up studies of larger source samples are 
needed to test the statistical significance of the hypothesis by 
SSHIO. 

5.2. Depletion and deuteration 

The CO depletion factors derived in this work are small, only in 
the range ~ 0.6 - 2.7. Thus, CO molecules do not appear to be 
significantly depleted, if at all, in the studied clumps. This con- 
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forms to the fact that the clumps show signs of star formation ac- 
tivity, such as outflows, which presumably release CO from the 
dust grains back into the gas phase (Appendix B). The presence 
of outflow-shocked gas is revealed by the SiO line wings (par- 
ticularly towards 118151 MM2) and very high SiO abundances 
found in the present work. Also, the temperature in the envelopes 
of 118102 MMl and 118151 MM2 is, likely due to heating by 
embedded YSOs, slightly higher than the CO sublimation tem- 
perature of 20 K (e.g.. lAikawa et al.~2 008 ). Moreover, relatively 
high cosmic -ray ionisation rates towards the clumps were found. 
Thus, the cosmic-ray impulsive heating could also play a role in 
returning CO into the gas phase (Has egawa & Herbst 1993| ). 

However, the derived CO depletion factors represent only 
the average values along the line of sight within the 27'.'8 beam. 
Because the studied clumps lie at kiloparsec distances from us, 
the observations could be "contaminated" by non-depleted gas 
along the line of sight dFontani et al. 2006 ). It is also uncertain 
whether the "canonical" CO abundance used to calculate /d is 
exactly correct. Its value is known to vary by a factor of ~ 2 bet- 
ween diff'erent star-forming regions (Lacy et al. 1994V Finally, 
the H2 column densities derived from dust emission suffer from 
the uncertainty of dust opacity, which may be a factor of 2-3. 
These uncertainties are likely to explain the unsual values of 
/d<1. 

For comparison, Fontani et al. (2006) determined the amount 
of CO depletion in HMYSOs. They found CO depletion factors 
in the range 0.4-35.8 with a median value of 3.2. Thomas & 
Fuller (2008) found depletion factors < 10 for their sample of 
HMYSOs and deduced that the sources are young, a few times 
lO^* yr. 

At high angular resolution, it is possible to identify more CO 
depleted regions within the clumps. Zhang et al. (2009) studied 
two massive clumps in the IRDC G28. 34-1-0. 06 at ~ 1" reso- 
lution with the SMA. They found CO depletion factors up to 
~ 10^ - lO"* towards the two clumps, where the largest values 
were found in the cores within one of the clumps. Similarly, 
with the ~ 3" resolution PdBl observations, Beuther & Henning 
(2009) found a large CS depletion factor of ~ 100 in IRDC 
19175-4. 

The degree of deuteration in our clumps are ~ 0.0002-0.014 
in HCO+ and ~ 0.002 - 0.028 in N2H+. Compared to the ave- 
rage cosmic D/H ratio of ~ 1.5 x 10"^ (Sect. 4.4), the derived 
values are ~ 10' - lO^* times larger Relatively low deuteration 
degrees in some of our sources are consistent with the oberved 
low CO depletion factors. CO is the main destroyer of H3 and 
HaD^, and thus its presence can lower the deuterium fractio- 
nation (e.g.. lCaselli et al. 20081 1. For example. Swift (2009) esti- 
mated the fractional abundance of ortho-H2D^ in two IRDCs of 
only ~ 3 - 5 X 10 which is up to three orders of magnitude 
lower than found in low-mass cores (e.g., ICaselh et al. 20081 
IFriesen et al. 2010 1). From a theoretical point of view, if there are 
no depletion gradients in the source, one would expect to see the 
equality /?d(HCO+) = /?d(N2H+) ( jRodgers & Charnley 2001bj . 
The fact that we found somewhat lower values for iJD(HCO^) 
than /?d(N2H^) could be related to differential depletion of 
molecules, and radial density gradients (' Caselli et al. 200 2 ). For 
high-mass star-forming clumps, the density profile is found to be 



of the form n(r) oc r 



-(1.6±0.5) 



(e.g., Beuthe r et al. 2002a). CO is 



the parent species of HCO^, but it destroys the N2H^ molecules. 
Deuteration proceeds most efficiently in regions where CO is 
mostly frozen onto dust grains. In the warmer envelope layers 
where CO is not depleted, HCO^ can have a relatively high 
abundance, i.e., a lower value of ^d(HCO^) compared to that 
of N2H^ ( Emprechtinger et al. 2009). In general, the amount of 



molecular deuteration is expected to decrease during protostel- 
lar evolution because of internal heating of the surrounding en- 
velope (e.g., Emprechtinger et al. 2009]). Indeed, there is a de- 
creasing trend in 7?d(HCO^) with the gas kinetic temperature as 
shown in Fig. |9] When the temperature becomes > 20 K, reac- 
tion 2 in Table [8]can proceed in both directions, and thus the de- 
struction rate of H2D^ increases. Consequently, other deuterated 
molecules that form via reactions with H2D^ start to decrease in 
abundance. 

For comparison, Fontani et al. (2006) determined 7?d(N2H^) 
in HMYSOs and found the values in the range < 0.004 - 0.02, 
with an average value of ~ 0.015. Roberts & Millar (2007) de- 
rived the upper limits to /;d(N2H+) of < 0.002 and < 0.006 
for the hot cores G34.26 and G75.78, respectively. Our source 
118102 MMl appears similar in this regard. Chen et al. (2010) 
found that 7?d(N2H+) = 0.017 - 0.052 in three cores within 
the IRDC G28.34h-0.06. Moreover, fliey found that 7;d(N2H+) 
is lower in the more evolved stages of protostellar evolution. 
Recently, Fontani et al. (201 1) found the values of /;d(N2H+) = 
0.012 - 0.7, 0.017-< 0.4, and 0.017-0.08 for their sample of 
high-mass starless cores, HMYSOs, and UC Hu regions, re- 
spectively. The above results are quite similar to our findings. 
Another study of deuteration in IRDCs is that by Pillai et al. 
(2007) who derived NH2D/NH3 column density ratios in the 
range 0.005-0.386 for their sample of IRDC clumps. More re- 
cently, PiUai et al. (201 1) found [NH2D/NH3] ratios in the range 
0.06-0.37 towards IRDCs. The [NH2D/NH3] ratios obtained by 
Pillai et al. are even higher than those observed in low-mass 
dense cores (see lPillai et al. 20Tn and references therein). 

Again, high resolution interferometric observations could 
reveal substructures within the clumps with higher degree of 
deuteration. For example, Fontani et al. (2006) found an ave- 
rage value /?d(N2H^) ~ 0.01 towards the high-mass star- 
forming region IRAS 05345-1-3157, but at high angular reso- 
lution (~ 3") they resolved it into two N2D^ condensations 
each with ten times higher deuteration degree [/?d(N2H^) = 
0.11; IFontani et al. 2008J . Such high values are comparable 



to those found in low-mass starless cores | Crapsi et al. 2005 
j;D(N2H+) < 0.02 - 0.44] and Class protostellar envelopes 



JEmprechtinger et al. 2009[ 7?d(N2H+) < 0.029 - 0.271]. For 
their sample of low-mass dense cores, Caselli et al. (1998) found 
/?d(HCO'^) to lie in the range 0.025-0.07. For a sample of more 
massive cores, Bergin et al. (1999) found that most sources have 
values of /?d(HCO^) < 0.035. This trend conforms to our re- 
sults towards more massive clumps. It is possible that the cold 
phase during which the deuterium fractionation takes place is so 
short for HMYSOs that very high deuteration degrees are not 
reached (Rob erts & Millar 2 007). On the other hand, the high 
[NH2D/NH3] ratios found by Pillai et al. (2007, 201 1) could in- 
dicate a different production mechanism of deuterated ammonia 
compared to those of N2H^ and HCO^ which are both purely 
produced in the gas phase. 

One caveat to the deuteration analysis is that the FWHM 
linewidths of the detected N2H^ transitions are significantly 
larger than those of N2D+ [Table g] Col. (4)]. Similarly, the 
linewidths of Hi3CO+(3 - 2) and DCO+(3 - 2) are quite dif- 
ferent compared to each other towards some of the sources. The 
above transition pairs have similar critical densities, and the lines 
were observed at comparable spatial resolutions. Therefore, the 
difference between the linewidths suggests that the correspond- 
ing transitions may not be tracing the same gas component. In 
this picture, the calculated deuteration degrees should only be 
interpreted as an average values along the line of sight. 



13 



Miettinen et al.: Massive clumps in IRDCs 



0.01 5 



8 0.010 



0.005 



0.000 



14 16 18 20 
T. [Kl 



24 



Fig. 9. A^(DCO^)/A^(HCO^) ratio versus gas kinetic tempera- 
ture. As Tidn increases, the degree of deuterium fractionation de- 
creases. 



5.3. lonisation 

We have calculated the average ionisation degree along the line 
of sight towards the clumps by utilising three different methods. 
The summed abundances of the ionic species yield lower lim- 
its to x(e) in the range ~ 0.7 - 12.3 x 10"**. A very simple 
chemical scheme including HCO^ and N2H^ yielded another 
estimates for the lower limits, about 0.3 - 5.6 x 10"^. Thirdly, 
we used the 7?d(HCO^) values to derive an upper limit to frac- 
tional ionisation, and found the range x(e) - 0.2 - 29.2 x 10"^. 
To our knowledge, these are the first reported estimates of the 
ionisation degree in IRDCs made so far As discussed below, the 
derived lower limits are likely to be more realistic estimates of 
x{e), whereas the upper limits appear to be very high. The x(e) 
values should be taken as rough estimates because of the sim- 
plicity of the analytical model we have utilised. For example, 
our analysis is based on chemical equilibrium which may not be 
vaUd Caselli etal. 19981 [Lintott & Rawlings 2 006 ). Other fac- 
tors that cause uncertainties are the possible variation of CO de- 
pletion and fractional molecular abundances along the line of 
sight. Electron abundance depends strongly on the depletion of 



CO because CO destroys ionic species such as and N2H^. In 
this regard, one could think that the lower limit to x(e) derived 
by summing up the abundances of ionic species is the most reli- 
able one. However, there is also a possibility that the line emis- 
sion from different molecules originate in different layers of the 
clump. Indeed, the different linewidths between N2H^(3 -2) and 
N2D^(3 - 2) point towards such a possibility (Sect. 5.2). In this 
case, the physical meaning of the summed abundance of ionic 
species is unclear Furthermore, the fact that we have no infor- 
mation about the abundances of some important molecular ions, 
such as Hl^ and H30^, nor the metal ions, hampers the determi- 
nation of .Ji:(e). Metal ions play an important role in the ionisa- 
tion balance in the outer layers of dense clouds (Caselli 20021 ). 
The fact that the sources show low degrees of deuteration con- 
forms to the fairly high ionisation levels of gas. However, de- 
tailed chemical modelling is needed to perform a more accurate 
analysis of the ionisation degree. 

The cosmic -ray ionisation rates of H2 we found are in the 
range ~ 1 x 10"'^ - 1 x 10"'^ s"'. Six of the seven sources 
have estimated values between ~ 1 x 10"'^ and ~ 5 x 10""' s"', 
but the highest rate of ~ 10"'^ s"' is seen towai'ds 118102 MMl. 
It should be emphasised that the values were derived from 
the lower limits to x{e). However, as discussed above, we be- 
lieve that these estimates are more accurate than those utilising 



the upper limits to x(e). Also, the appUcability of the volume- 
averaged H2 number densities used to compute is uncertain 
because all the detected molecular-ion transitions have a much 
higher critical density (Hez areh et al. 20()8b . The standard rela- 
tion between the electron abundance and the H2 number density 
is x(e) ~ 1.3 X 10"^n(H2)-^^2 dMcKee 19891 ). This is based on 
the pure cosmic -ray ionisation with the rate 1.3 x 10"'^ s"' and 
includes no depletion of heavy elements. By utilising the den- 
sities derived in the present work, 1.1 - 18.5 x 10"^ cm"-*, the 
standard relation yields the values x(e) ~ 3 - 12 x 10"^. These 
are comparable to the estimated lower limits to x(e). 

How do our results compare to those found in other sources ? 
Caselli et al. (1998) determined x(e) in a sample of 24 low-mass 
cores consisting of both starless and protostellar objects. Their 
analysis was based on observations of CO, HCO^, and DCO^, 
and the resulting values were in the range 10"^ to lO"*". They ar- 
gued that the variation in x(e) among the sources is due to vari- 
ations in metal abundance and Williams et al. (1998) used 
observations of C'^O, H'^CO+, and DCO+ to determine the val- 
ues lO"'''^ < x(e) < lO"'''^ in a similar sample of low-mass cores 
as Caselli et al. (1998), but they used a slightly different analysis. 
Applying the same analysis as Williams et al. (1998), Bergin et 
al. (1999) found the ionisation levels of 10"^^ < x(e) < 10"*' '' 
towards more massive cores (in Orion) than to those studied by 
Williams et al. The most massive sources were found to have the 
lowest electron abundances, .Ji:(e) < 10"^. This conforms to the 
results by de Boisanger et al. (1996) who found that .Ji:(e) ~ 10"^ 
in the massive star-forming regions NGC 2264 IRSl and W3 
IRS5. Hezareh et al. (2008) studied the high-mass star-forming 
region DR21(OH) in Cygnus X, and found that x(e) = 3.2x lO l 
These are consistent with our lower limits to .Ji:(e) in massive 
clumps within IRDCs. The very high upper limits to x(e) we 
have obtained, particulai'ly that in 118102 MMl (~ 10""*), re- 
semble those found in PDRs, where carbon can provide most of 
the charge and x(e) ^ x(C^) ( iGoicoechea et al. 2009i ). 

Observational results presented in the literature suggest re- 
latively large environmental variations in ^hJ^. van der Tak & 
van Dishoeck (2000) used H^ and H'^CO^ observations to con- 
strain the ionisation rate towards HMYSO envelopes, and ob- 
tained the best-fit value (n. = (2.6 + 1.8) x 10"'^ s"'. We 
found comparable rates towards three of our sources. Caselli 
et al. (1998) inferred the values of spanning a range of 
two orders of magnitude between 10"'* - 10"'^ s"' in low- 
mass cores. Some of this variation could be due to different 
cosmic -ray flux in the source regions. WUiams et al. (1998) de- 
duced a mean value of ^h, = 5 x 10"'^ s"'. On the other hand, 
H^ observations towards the diffure ISM suggest higher values 
in die range fa, ~ 0.5 - 1.2 x 10"'^ s"' jMcCaU et al. 20031 
llndriolo et al. 20071) . Also, towards the Galactic Centre clouds, 
the ionisation rate is found to be very high, fa, > 10"'^ s~' up 
to ~ IQ-'^ S-' (e.g., lOka et al. 20051 lYusef-Zadeh et al. 20071 
IGoto et al. 20081 ). However, such high rates can explain the high 
temperatures of ~ 70 - 100 K in the Galactic Centre clouds, 
whereas the dense YSO envelopes are typically much colder van 
der Tak et al. (2006) estimated that fa, ~ 4 x 10"'* s"' in the Sgr 
B2 envelope, whereas Hezareh et al. (2008) found a lower value 
of fa, = 3.1 X 10""* s"' towards DR21(OH). Observations of the 
Horsehead Nebula by Goicoechea et al. (2009) could only be 
reproduced with fa, = 7.7 + 4.6 x 10"'^ s"'. Despite aU the 
uncertainty factors in fa,, it seems possible that some of our 
clumps could be exposed to somewhat higher cosmic-ray flux. 



In the following, we have used the relation 2.3fH = l-5<rH2 to con- 
vert the ionisation rate per H atom, (h, into ■ 
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but the observed clump temperatures, ^ 20 K, suggest that is 
~ 10"'^ - 10 (see|Bergin et al. 1999 ). 



Table 10. The field-neutral coupling parameter (W), MHD cutoff 
wavelength (/Imin), and ambipolar diffusion timescale (tad)- 



Because (most of) our sources are associated with high-mass 
star formation, which could increase the local UV radiation field, 
the role played by photoionisation should also be considered. 
Photoionisation is expected to become negligible when the ex- 
tinction is Ay > 4 mag (McKee 1989). In units of A^(H2), thi s 
corresponds to about > 4.9 x 10^' cm"^ ( Vuong et al. 2003 ). 



Therefore, the ionisation in the studied clumps is very likely 
dominated by cosmic ray particles and the ambient UV radi- 
ation field is not important in this regard. If the clumps are 
inhomogeneous, however, external UV radiation could pene- 
trate more efficiently into them and enhance the ionisation level 
JBoisse 19901 ). Our sources also show high HCO^ abundances 
(~ 5 X 10"' - 9 X 10"^), whereas in the presence of signifi- 
cant UV radiation one would expect to observe lower values 
of x(HCO^) due to photoproduced electrons raising the HCO^ 
recombination rate ( IJansen et al. 1995) . On the other hand, the 
HCO^ abundance could be enhanced by the presence of proto- 
stellar outflows where shock-generated UV radiation is present 
(Rawlings et al. 2000). This could explain the high HCO^ abun- 
dance in 118102 MMl (Sect. B.l). 



5.4. Coupling between the ions and neutrals 

The knowledge of the gas ionisation degree enables to estimate 
the level of coupling between the neutral gas component and 
the pervasive magnetic field. The coupling strength can be quan- 
tified through the wave-coupling number, which is defined as 
the ratio of the clump radius to the MHD cutoff wavelength, 
W = R/A^in ( Myers & Khersonsky 1995) . The wavelength 
represents the minimum wavelength for propagation of MHD 
waves, and is given by A^in = jrvATin, where Va is the Alfven 
speed, and Tin is the ion-neutral collision time. When A < /tmin, 
the frequency of the wave is higher than the collision frequency 
between the ions and neutrals (v » vin). In this case, the neutral 
fluid component is decoupled from the magnetic-field dynamics. 
If the clump is much larger than /imn, i.e., W » 1, the magnetic 
field is strongly coupled to the neutral component, and MHD 
waves can stabilise the clump. 

We calculated the W number by utilising Eq. (6) in Williams 
et al. (1998): 



W 



<e)n(H2)i/2<(ri„v) 



5(1 +P) \ 
AnnmnG - IScrlj [n{V{2)R^] ] 



1/2 



(14) 

where (cTjnv) is the ion-neutral collision rate, p e [0, 1] is a 
parameter describing the level of turbulence (0 for a minimum 
amount of turbulent motions, 1 for a maximum level of turbu- 
lence), ji is the mean molecular weight per particle (taken to be 
2.33), G is the gravitational constant, and cr-j- = ^^jk^f^^jjimii 
is the one-dimensional thermal velocity dispersion. Following 
WilHams et al. (1998), we take <cri„v) = 1.5 x 10"'* cm^ s"'. We 
used as x(e) the summed abundance of ionic species, and thus 
the values of W should be taken as lower limits. We assumed that 
the clumps are turbulent (p - \), and clump radii were taken to 
be the effective radii listed in Col. (7) of Table |5] Note that the 
latter were also used to derive the H2 number densities. The com- 
puted values of W and A^^m are given in Table [10] The derived 
W values are in the range ~ 2-19, whereas the minimum wave- 
lengths lie between 0.018 and 0.174pc. Note that the clump radii 
range from 0.15 to 0.34 pc, and are mostly comparable to each 
other Thus, the large variation of the W values cannot (solely) 



Source 


W 


^min 


Tad" 






[pc] 


[10' yr] 


IRDC 18102-1800 MMl 


18.5 


0.018 


37/162 


G015.05+00.07 MMl 


2.1 


0.072 


3/12 


IRDC 18151-1208 MM2 


2.6 


0.065 


2/9 


G015.31-00.16MM3 


2.2 


0.174 


7/30 


IRDC 18182-1433 MM2 


6.0 


0.044 


16/65 


IRDC 18223-1243 MM3 


4.9 


0.037 


5/24 


ISOSS J18364-0221 SMMl 


8.9 


0.036 


12/54 



Notes. The first Tad value was estimated from the computed value 
of W, and the second one was calculated by assuming that HCO^ is the 
dominant ionic species (see text). 



be explained by different source sizes. Because W represents the 
ratio of maximum to minimum wavelength, the maximum wave- 
length corresponds to the clump radius. For comparison, Bergin 
et al. (1999) deduced that massive cores in Orion (comparable 
in size to our clumps) are characterised by the value W - 20. 
The MHD wave power transmission at the lowest frequency of 
wave propagation (i.e., the longest wavelength) rapidly drops 
below unity when W < 100 ( [Myers & Lazariarri998 ; their 
Fig. 1). Therefore, the values W < 100 imply that the coupling 
is "marginal", and that MHD waves have a quite narrow band of 
wavelengths to propagate above A^i„ . However, for 1 1 8 1 02 MM 1 
and J18364 SMMl the allowed wavelengths of MHD waves 
could extend from the effective radii ~ 0.2 - 0.3 pc down to 
Amin ~ 0.02 - 0.04 pc. This is expected to lead to a clumpy struc- 
ture of a medium (Myers & Khersonsky 1995} , just like is found 
to be the case in J 18364 SMMl which is fragmented into two 
subcores (see Sect. B.6). 

Another approach to investigate the coupling of the neutral 
gas and the magnetic field is to determine the ratio between 
the ambipolar diffusion (AD) timescale and free-fall timescale, 
T"AD/T"ff. In the process of AD, neutrals drift quasi-statically (un- 
der their own self-gravity) relative to ions and magnetic field 
( Mestel & Spitzer 1956| ). It can be shown, t hat the ratio TAp/Tg 
is approximately given by W ~ Tft /Ti,, (see lMouschovias 19911 
[Williams et al. 1998b . The free-fall timescales of our clumps. 



Tff = yl37T/32G{p) ~ 8 X 10** - 3 X 10^ yr, imply AD timescales 
in the range tad ~WXTff~2xlO^-4xl0^yr [see the first val- 
ues in Col. (4) of Table [To). We note that tad depends strongly 
on the ionic composition, and can be even much longer than the 
above rough estimates. Assuming that HCO^ is the dominant 
ionic species, the AD timescale is given by tad ~ 1.3 x 10'^Ji:(e) 
yr [see Eq. (5) in Walmsley et al. (2004)]. Using again the lower 
limits to .5t(e), we obtain timescales in the range Tad ~ 1-16x10^ 
yr which are about four times longer than the above estimates 
[see the second values in Col. (4) of Table [TOl . Therefore, on 
the scale probed by our single-dish observations, AD is not ex- 
pected to be a main driver of clump evolution unless it occurs 
on timescales much larger than 10^ yr. However, the situation 
may be different at smaller scales, where density is higher, and 
therefore x(e) is lower and tad is shorter 

6. Summary and conclusions 

We have carried out a molecular-line study of seven massive 
clumps associated with IRDCs. We used APEX to observe tran- 
sitions of C'^O, H'^CO^, DCO+, N2H+, and N2D+. The princi- 
pal aim of this study was to investigate the depletion, deuterium 
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fractionation, and tiie degree of ionisation in the sources. Our 
main results are summarised as follows: 

1. The CO molecules do not appear to be significantly de- 
pleted, if at all, in the observed sources. The largest CO deple- 
tion factor, ~ 2.7, is found towards G015.3 1 MM3 which is dark 
at MIR. In many sources, CO is Ukely to be released from dust 
grains into the gas phase due to shock evaporation. 

2. The deuteration degree in HCO^ was found to range from 
0.0002 to 0.014, whereas that in N2H+ lies in the range 0.002- 
0.028. These are smaller values than found in low-mass star- 
less cores and protostellar envelopes, but still significantly larger 
than the cosmic D/H-ratio ~ 10"^. The degree of deuteration in 
HCO^ appears to decrease with increasing gas temperature, as 
expected theoretically. This is hkely to reflect the evolutionary 
stage of the clump: at early stages the source is cold and as the 
evolutionary stage progresses the temperature rises. In the course 
of evolution the amount of deuterium fractionation decreases. 

3. For the first time, we have estimated the level of ionisa- 
tion in IRDCs. The ionisation degree is difiicult to determine ac- 
curately, especially from the analytical basis like in the present 
work. However, a lower limit to electron abundance appears to 
be x(e) > 10"** - 10"^, which is comparable to those seen in 
other star-forming regions. The cosmic-ray ionisation rate of H2 
were estimated to vary from ~ 10"'^ s~^ up to ~ 10"^^ s"^ in one 
target position. Similar estimates have been found in a variety of 
Galactic environments. 

4. Additional molecular species were detected towards some 
of the sources. In particular, ortho-c-C3ll2 was detected in four 
sources, and SiO towards three sources. This is the first reported 
observation of C-C3H2 towards IRDCs. The SiO lines showed 
large hnewidths indicative of outflows. Also for the first time 
towards IRDCs, the deuterated ethynyl, C2D, was detected in 
118151-1208 MM2. We estimate flie C2D/C2H column density 
ratio to be about 0.07, in good agreement with the results from 
other sources (Appendix A). 

5. The ambipolar-diffusion timescale imphed by the lower 
Umits to x(e) is typically ~several Myr Therefore, on the scale 
probed by our observations, clump evolution cannot be mainly 
driven by ambipolar diffusion unless it occurs on timescales 
much larger than 10^ yr. 

The present single-dish study cannot follow the chemical 
structure of the sources down to small scales where the evolu- 
tion of the massive young stellar object(s) is taking place (e.g., 
the case of J 18364 SMMl). Follow-up high-resolution studies 
are needed i) to resolve whether the clumps are fragmented into 
smaller subunits (as is known to be the case in J18364 SMMl); 
ii) to examine the kinematics (infall, outflow, relative motions) 
of the subfragments, and, importantly, their chemical partition 
(such as depletion and deuteration); and Hi) to locate the detected 
complex molecules (envelope or hot core ?). In this regard, 
follow-up observations with ALMA seem particularly worth- 
while. Finally, detailed chemical models are needed to better un- 
derstand the observed molecular abundances, and how these are 
related to the timescale of the high-mass star formation process. 
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Fig. 1. Spitzer/IRAC 8-jum images overlaid with contours of (sub)mm dust continuum emission (SCUBA 850 fim, Bolocam 1 . 1 mm, 
MAMBO 1.2 mm). The colour bars indicate the 8-/im intensity scales in units of MJy sr"'. The contour step-size is 0.1 Jy beam"' 
in all cases except G015.31 MM3 where it is 0.02 Jy beam"'. The contours go as follows: from 0.02 to 0.14 Jy beam"' for G015.31 
MM3; from 0.1 to 0.7 Jy beam"' for G015.05 MMl and 118151 MM2; from 0.2 to 0.7 Jy beam"' for J18364 SMMl; from 0.2 to 
0.9 Jy beam"' for 118182 MM2; from 0.2 to 2.3 Jy beam"' for 118102 MMl; and from 0.3 to 0.7 Jy beam"' for 118223 MM3. The 
crosses indicate the positions of our molecular-line observations. The box symbol towards 118102 MMl indicates the position of the 
6.7 GHz Class II CH3OH maser dBeuther et al. 2002bl l. For 118151 MM2, the box symbol shows the position of the 84.5 GHz Class 
I CH3OH maser (Marseille et al. 2010bl l: it coincides with our target position. The circles indicate the positions of the 22.2 GHz 
H2O masers associated with 118151 MM2 (Be uther et al. 2002bl l and GO 15. 05 MMl (|Wang et al. 2006| . The mm peak positions of 
the northern and southern cores within the clump J18364 SMMl are indicated jHennemann et al. 20091 1 . The scale bar in each panel 
corresponds to 0.5 pc. 
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Fig. 2. Smoothed C"0(2 - 1), H^^CO+CS - 2), DCO+(3 - 2), N2H+(3 - 2), and N2D+(3 - 2) spectra (fop to bottom) towards the 
clumps (left to right). Overlaid on the spectra are the hf-structure fits. The relative velocities of each individual hf components in 
each observed transition are labelled with a short bar on the spectra towards 118102 MMl; these are also shown on the and 
N2H+ spectra towards G015.31 MM3. The red vertical fine on the N2D+ spectrum towards G015.31 MM3 shows the Vlsr - 30.81 
km s"' of the clump as measured from N2H^(1 - 0) by SSK08. Note that N2H^(3 - 2) observations were carried out only towards 
three sources and thatN2D+(3 - 2) was not detected towards G015.31 MM3, 118182 MM2, and J18364 SMMl. 
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Fig. 4. orf/!o-c-C3H2(33.o - 22,i) spectra overlaid with Gaussian fits. 
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Fig. 5. SiO(6 - 5) spectra overlaid with Gaussian fits. The vertical line in each panel indicates the Vlsr of the clump as measured 
from N2H+(1 - 0) by SSK08, i.e., 22.36, 29.74, and 45.73 km s"' for 118102 MMl, 118151 MM2, and 118223 MM3, respectively. 
The velocity range in each panel is wider relative to that of Fig.|2]for a better illustration. 
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Fig. 6. C"^0(2 - 1) lines at 219 560.3568 MHz (JPL) arising from the image sideband. 
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Fig. 8. The range of molecular column densities (upper panel) and fractional abundances (lower panel) for our sources as com- 
pared with similar samples of other studies. The y-axes are shown on a logarithmic scale. The abbreviation TW refers to 'this 
work'. The references are coded as follows: FCC06=Fontani et al. (2006), SSH10=Sakai et al. (2010), BLL95=Butner et al. 
(1995), GPB09=Gibson et al. 2009, SSK08=Sakai et al. (2008), BH09=Beuther & Henning (2009), CLSlO^Chen et al. (2010), 
FPCll^Fontani et al. (2011), RBP06=Ragan et al. (2006), TF08=Thomas & Fuller (2008), VLHl l^Vasyunina et al. (2011). The 
'^CS abundance of ~ 6.2 x 10"'' we derived towards 118102 MMl is shown by a slightly stretch stock for a better illustration. 
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Appendix A: Other detected transitions 

Besides the molecular species discussed in Sect. 5.1, the sources 
118102 MMl, 118151 MM2, 118182 MM2, and 118223 MM3 
show emission lines also from other species, some of which are 
complex organics. These include 0-C-C3H2, C2D, C'^O (seen 
in the image band), and possibly CNCHO (in the image band), 
CH3NH2 (blended with OCS), CHiCHC'^N, and CH3COCH3. 
In this Appendix, we discuss each of these species (except C'^O) 
and their derived column densities and abundances. 

Cyclopropenylidene (c-C-iHi)- We have detected the ortho 
form of cyclic-C3H2 in four sources, and found the column den- 
sities and abundances in the range ~ 7.0 x 10'^ -7.7 x lO'"* cm"^ 
and ~ 1 .9x 10 1.2x 10 ^ respectively. The C-C3H2 molecule 
has been detected in several different Galactic sources but to our 
knowledge, this is the first reported detection of this molecule to- 
wards massive clumps associated with IRDCs. For comparison, 
the beam-averaged C-C3H2 column densities and abundances of 
~ 4 - 9 x 10'^ cm~^ and 3 - 7 x 10"'^ have been found in Sgr 
B2 dTurner 19911 INummeli n et al. 20 00). These column densi- 
ties are comparable to those in our clumps, but the abundances 
appear lower in Sgr B2. 

The major formation channel of C-C3H2 is €3!!^ -n H2 ^ 
c - C3H+; c - C3Ht -H e- ^ c - C3H2 (e.g., IPark et al. 2006b . 
The very low rotational excitation temperatures we have derived, 
3. 1^.8 K, suggest that C-C3H2 emission mainly comes from the 
cool and relatively low-density envelope dTurner 19911 1. Also, 
higher electron abundance in the outer layers can promote the 
formation of C-C3H2. This is consistent with the derived high 
fractional abundances. 

Cyanoformaldehyde or formyl cyanide (CNCHO). The 
CNCHO transition is possibly seen in the image sideband to- 
wards 118102 MMl. Because of attenuated intensity, we cannot 
determine the column density of this molecule. Remijan et al. 
(2008) determined the values A?(CNCHO) =1-17x10'"* cm^^ 
and jc(CNCHO) = 0.7 - 1 1 x lO '' towards Sgr B2(N). They sug- 
gested that CNCHO is likely formed in a neutral-radical reaction 
of formaldehyde (H2CO) and the cyanide (CN) radical. 

Methylamine ( CHt,NH2)- Rotational levels of CH3NH2 un- 
dergo the so-called A- and /s-type transitions due to the inter- 
nal rotation of the CH3 -group. We have a possible detection 
of the ii-type CH3NH2 transition, and we derived the values 
A^(CH3NH2) = 5x10''' cm-2 andx(CH3NH2) = 5.1 +0.6 x 10^1 
However, our CH3NH2 line is blended with OCS line, and thus 
these values should be taken with caution. Beam-averaged co- 
lumn densities of A^(C H3NH2) ~ 1 X 10'"* - 8 X 10'^ cm"^ have 
been found in Sgr B2 ( Turner 199 11 INununehn et al. 20001) . The 
origin of CH3NH2, and nitrogen-bearing organics in general, 
could be in the evaporation of ice mantles, indicative of a hot- 
core chemistry (e.g. Rodgers & Charnley 2001a 1. 

Carbonyl sulfide (OCS). The high-7 transition of OCS we 
have possibly detected (blended with CH3NH2) implies a co- 
lumn density and fractional abundance of 2.3 + 0.8 x 10'^ cm"^ 
and 2.4 ± 0.9 x 10^^, respectively. For example, Qin et al. (2010) 
detected OCS in the high-mass star-forming region G19. 61-0. 23. 
Assuming Trot = E^jk^, they derived the beam-averaged column 
density and fractional abundance of 2.2 ± 0.1 x 10"' cm"^ and 
2.7 + 0.1 X 10"**, respectively. Note that Qin et al. achieved a 
much higher angular resolution with their SMA observations, 
and thus their large column density value could partly be due to 
filtering out the extended envelope. Our coarser angular resolu- 
tion probably causes a significant beam dilution. SSHIO derived 
OCS column densities of 1.2 - 5.5 x lO'** cm"^ for their sample 



of 20 massive clumps associated with IRDCs, i.e., lower than 
found here towards 118102 MMl. 

The OCS molecules form on grain surfaces through the ad- 
dition of S atom to CO, or via the O atom addition to CS (e.g.. 



Charnley et al. 2004 1. Sohd-state OCS has been detected to- 



wards high-mass star-forming regions by, e.g, Gibb et al. (2004) 
SSHIO suggested that OCS is released from the grains into the 
gas phase through protostellar shocks. This could well be the 
case in 118102 MMl (Sect. B.l). 

Deiiterated ethynyl (C2D). We have made the first detec- 
tion of C2D towards IRDCs. The column density and abundance 
estimated from the N - 3 - 2 transition are ~ 1.5 x lO'-' cm"^ 
and 4.2 + 0.5 x 10"'°, respectively. SSHIO detected the normal 
isotopologue C2H(A^ =1-0) towards 118151 MM2. The C2H 
column density they derived, ~ 2.2 x lO'** cm"^, together with 
A^(C2D) derived by us, suggest a deuteration degree of ~ 0.07 
in C2H. Vrtilek et al. (1985) detected the CiDiN = 2-1) 
transition near the Orion-KL position, and obtained the co- 
lumn density ~ 1.8 x 10'^ cm"^. Moreover, they derived the 
A^(C2D)/A^(C2H) ratio of 0.05. These are very similar to what 
we have found. Quite similarly, Parise et al. (2009) found upper 
Umits of A^(C2D) < 2.5 x lO'^ cm'^ and x(C2D) < 2 x 10"'" in 
a clump associated with the Orion Bar. 

The formation of C2D is believed to take place in the gas 
phase through the route CH3 — > CH2D+ C2D (see, e.g., 
IParise et al. 2009] l. The A?(C2D)/A^(C2H) ratio we have obtained 
is comparable with those predicted by the low-metal abundance 
model by Roueff et al. (2007) at temperature around 30^0 K. 
We note that the C2H abundance, ~ 6 x 10 calculated from the 
observed deuteration degree (0.07) is comparable to the values 
x(C2H) = 2.5 X 10"^ - 5.3 x 10"** recently found by Vasyunina 
et al. (201 1) towai-ds IRDCs. 

Acrylonitrile or vinyl cyanide, ^^N isotopologue 
(CH2CHC^^N). The Weeds modelUng suggests that there 
is a CH2CHC'^N transition blended with the hf group of the 
detected C2D line. Thus, reliable column density estimate can- 
not be performed. Assuming that the detected line is completely 
due to CH2CHC'^N emission, we derive very high values of 
2.5 X 10'* cm"^ and ~ 7 x 10"^ for the column density and 
fractional abundance. For comparison, a column-density upper 
limit of 3 X 10'^ cm"^ towards Sgr B2(N) was derived by Miiller 
et al. (2008). 

If present, this molecule would indicate a hot-core che- 
mistry. The main isotopologue C2H3CN is expected to form 
through gas-phase reactions after the ethyl cyanide (C2H5CN), 
forming on dust grains, evaporates into the gas phase 
dCaselh et al."7993T l. 

Acetone (CH3 COCH3). The first detection of CH3COCH3 
in the interstellar medium was made by Combes et al. (1987) 
towards Sgr B2. They found the column density and fractional 
abundance of this molecule to be 5 x lO'"* cm"^ and 5 x 10"". 
These are significantly lower than we have estimated towards 
118182 MM2 (2 x 10'^ cm^^ and ~ 4 x 10"^). Combes et 
al. (1987) found that the CH3COCH3 abundance is about 1/15 
of its precursor molecule CH3CHO (acetaldehyde), and sug- 
gested the formation route of acetone to be CH3 + CH3CHO — > 
(CH3)2CHO+ + hv; (CH3)2CHO+ -n e' ^ CH3COCH3 + H. 
However, Herbst et al. (1990) showed that this radiative as- 
sociation reaction is likely too slow to be consistent with the 
observed abundance in Sgr B2. The chemistry behind the for- 
mation of acetone is not clear. It could be due to some other 
gas-phase ion-molecule reactions, or due to grain chemistry 
dHerbst et al. 19901) . 
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Friedel et al. (2005) found that acetone in Orion BN/KL 
is concentrated towards the hot core. They derived the beam- 
averaged column densities of ~ 2 - 8 x 10"' cm"^. More recently, 
Goddi et al. (2009) found the column density A^(CH3COCH3) = 
5.5 X 10"' cm"^ in Orion BN/KL, in agreement with the Friedel 
et al. results. The very high column densities and abundances 
of CH3COCH3 found in Orion BN/KL and in 118182 MM2 in 
the present work require some other reaction pathway(s) than the 
above radiative association reaction to be efficient. Grain surface 
and hot-core gas-phase chemistry may both play critical roles 
dOarrod et al. 20081 1. 

The detection of complex molecules in clumps associated 
with IRDCs, indicating the presence of hot cores, supports the 
scenario that high-mass star formation can take place in these 
objects. For example, Rathborne et al. (2007, 2008) found that 
the clumps G024.33+00.1 1 MMl and G034.43+00.24 MMl 
are both likely to contain a hot molecular core. On the other 
hand, complex organics could also be ejected from grain man- 
tles through shocks, and cosmic rays heating up the dust can 
also have some effect (e.g., Requena- Torres et al. 2008 and re- 
ferences therein). 



Appendix B: Discussion on individual sources 

B.1. IRDC 18102-1800 MM1 

The clump 118102 MMl is the warmest (21.3 K) and most mas- 
sive (~ 360 Mq) source of our sample. It is associated with 
Spitzer point sources at 8 and 24 /vm, and high-mass star forma- 
tion is taking place within it as indicated by the presence of the 
6.7 GHz Class II CH3OH maser (Beuther et al. 2002b ). Among 
our sample, the lowest degree of deuteration in both HCO^ and 
N2H+ is found for 118102 MMl. On the other hand, the source 
shows the highest degree of ionisation. 

Fuller et al. (2005) detected central dips and red asymmetries 
in the spectral Unes HCO+(l - 0), HCO+(4 - 3), N2H+(1 - 0), 
and H2CO(2i 2 - li,i). A similar line profile is seen in the 
7 = 3-2 transition of N2H^ in the present study. These profiles 
indicate the presence of expanding or outflowing gas (see, e.g., 
IPark et al. 20001 and references therein). Beuther & Sridharan 
(2007) detected very broad (42.2 km s"' down to zero intensity) 
SiO(2 - 1) wings towards 118102 MMl, indicative of bipolar 
outflows. This source showed the second broadest SiO line in 
the sample of Beuther & Sridharan (2007). Also, the SiO(2 - 1) 
line detected by SSHIO towards 118102 MMl was very wide, 
Av = 13 ± 1 km s The SiO(6 - 5) line detected in the present 
study has Av ^ 4.7 km s ' , and the width down to zero intensity 
is ~ 10 km s"'. We have also detected the '^C isotopologue of 
CS in this source. SSHIO suggested that CS could originate in 
the shock evaporation of grain mantles or radiative heating. 

Beuther & Sridharan (2007) detected CH3CN(6jf - 5k) and 
CH30H(5;f - 4k) Hnes towards 118102 MMl, and derived the 
abundances of 3 x 10 " and 4 x 10 " for CH3CN and CH3OH, 
respectively. SSK08 detected CH30H(7^ -6/r) and HC3N(5 -4) 
lines in this source, and SSHIO detected C2H(A^ =1-0) and a 
hint of CH30H(2i - li) A . We note that SSK08 and SSHIO did 
not detected the lines of CCS(43 - 32), SO(22 - 1 1), or OCS(8 - 
7) in their surveys. The upper limit they derived for the OCS 
column density, < 2.7 x lO'"* cm"^, is about 8.5 times less than 
the A^(OCS) value we have obtained. The observational results 
cumulated so far indicate hot-core chemistry in 118102 MMl. 
This conforms to the fact that this clump is giving birth to a high- 
mass star(s), and possibly through disk accretion as indirectly 
suggested by the outflow signatures. The associated hot core is 



likely to be in its later stages of evolution because it is associated 
with a methanol maser (cf. lRathborne et al. 20081 1. 



B.2. G01 5.05+00.07 MMl and G015.31-00.16 MM3 

G015.05 MMl is the lowest mass (63 M©) and G015.31 MM3 
is the second coldest (13.7 K) clump of our sample. Both clumps 
are dark in the Spitzer 8 and 24 yum images. The highest deutera- 
tion degree in N2H+ (0.028) was found towards G015.05 MMl, 
whereas G015.31 MM3 shows the largest CO depletion fac- 
tor (2.7 ) in our sample. G 015.05 MMl is associated with H2O 
maser ( Wang et al. 2006[ l, indicative of star formation activity. 
Rathborne et al. (2010) derived the following properties for 
G015.05 MMl from the broadband SED: Tdust = 11 - 36 K, 
L = 15.5 - 362 Lq, and M = 35 - 158 Mq (scaled to the re- 
vised distance 2.6 kpc). The values Tkm = 17.2 K and M = 63 
M0 derived in the present study lie at the low end of the above 
values of temperature and mass [at the derived clump densities, 
it is expected that Tki,, - Tdust ( Goldsmith & Langer 1978l l]. 

SSK08 barely detected the HC3N(5 - 4), CH30H(7/f - 6;^), 
and CCS(43 - 32) lines towards G015.05 MMl and G015.31 
MM3. Actually, all the CH3OH detected objects in the survey 
by SSK08 are associated with the Spitzer lA-fim sources. 

Both GO 15. 05 MMl and GO 15. 31 MM3 are likely to be in a 
very early stage of evolution. Moreover, both sources are mas- 
sive enough to allow high-mass star formation. GO 15. 31 MM3 
could represent or host the so-called high-mass prestellar core, 
whereas GO 15. 05 MMl could be slightly more evolved with 
H2O maser emission but still lacking IR emission at 8 and 24 
fim. The small wave-coupling number of W ~ 2 for G015.31 
MM3 suggests that "magnetic turbulence" is not able to frag- 
ment the clump into smaller pieces, strengthening the possibility 
that it is a massive prestellar "core". 



B.3. IRDC 18151-1208 MM2 

The 118151 MM2 clump is dai'k in the Spitzer 8- 
pm image (24 pm not available). It is associated with 
H2O (Ee uther et al. 2002b) and Class I CH3OH masers 
(IMarseille et al. 2010b ). This clump has the highest volume- 
averaged H2 number density (~ 1.9 x 10^ cm"^) among our 
sources. 

Beuther & Sridharan (2007) detected ±s broadest SiO(2 - 1) 
wing emission (65 km s"' down to zero intensity) towards 
118151 MM2 in their sample. Also, the 7 = 2 - 1 and 3 - 2 SiO 
lines detected recently by Lopez-Sepulcre et al. (201 1) are very 
broad, i.e., FWZP= 84.3 and 103.1 km s"', respectively. The 
SiO(6 - 5) line we detected is also very broad with the FWHM 
46.8 km s"'. The outflow activity within the clump was con- 
firmed by Marseille et al. (2008) who, for the first time, found 
that 118151 MM2 is driving a CO outflow and hosts a mid-IR- 
quiet, possibly a Class 0-Uke HMYSO (cf. iMotte et al. 20071 ). 
Marseille et al. (2008) modelled the dust continuum emission 
(SED) of 1 1 8 1 5 1 MM2 and found that the bolometric luminosity, 
mass, and the mean temperature of the source are Lboi = 2190 
Lq, Mgas = 373 + 81 Mq, and (T) = 19.4 + 0.2 K (scaled to 
the revised distance 2.7 kpc). Within the errors, these mass and 
temperature values are comparable to the values derived in the 
present paper. Besides the Class I methanol maser tracing the 
molecular outflow, Marseille et al. (2010b) detected blue asym- 
metry in CH30H(5_i,5 - 4o,4) E, indicating infall motions. The 
22 GHz H2O maser is probably excited in the outflow shocked 
gas (cf. Furuya et al. 201 1 1. 
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Beuther & Sridharan (2007) also detected CH3CN(6;f - 5;^) 
and CH30H(5jf - 4k) lines towards 118151 MM2, which is a 
sign of hot-core chemistry, and they derived the abundances of 
8 X 10 " and 6 x 10 "' for CH^CN and CH3OH, respectively. 
SSK08 detected CHaOHC?;^ - 6k) lines towards 118151 MM2, 
but not CCS(43 - 32) or HC3N(5 -4). In their line survey, SSHIO 
detected C2H(A? =1-0), but not SO(22 - li), OCS(8 - 7), or 
CH30H(2i - 1 1) A" lines. The C2H column density they derived, 
~ 2.2 X 10'"* cm-^, together with the value A^(C2D) = 1.5 x 10'^ 
cm"^ derived by us, suggest a deuteration degree of ~ 0.07 in 
C2H. Marseille et al. (2008) concluded from their modelling of 
CS transitions that CS is depleted in 118151 MM2. We derived 
only a small CO depletion factor of ~ 1.6 for this clump, and 
low degrees of deuteration in HCO^ and N2H^, namely 0.3% 
and 1%, respectively. Also, this source shows the lowest lower 
limit to ionisation degree, only 7 x 10 

B.4. IRDC 18182-1433 MM2 

The filamentary clump 118182 MM2 is associated with Spitzer 
8 and 24-yum sources. This source shows the highest degree of 
deuteration in HCO+ (0.014). 

Beuther & Sridharan (2007) detected CH30H(5jf -4;^) lines 
towards 118182 MM2, and derived the CH3OH abundance of 
2.1 x 10-'°. SSK08 did not detect CCS(43 - 32), CH^OHilK - 
6k), or HC3N(5 - 4) lines towards this source in their survey. 
Our tentative detection of the O-bearing species CH3COCH3 
indicates hot-core chemistry within the clump. Moreover, O- 
bearing species are sign of the early stage of chemical evolution 
(e.g.. lShiao et al. 2010l l. This conforms to the presumable young 
age of the clump as it is associated with IRDC. Interestingly, 
the nearby clump 118182 MMl, which is associated with the 
HMYSOIRAS 1818 2-1433, also appears to contain a hot core 
dBeuther et al. 2006l l. 

as. IRDC 1 8223-1 243 MM3 

118223 MM3, part of a long filamentary IRDC, is a high-mass 
clump harbouring an embedded accreting low- to intermediate- 
mass protostar which could evolve to a high-mass star at 
some point in the future (see IBeuther & Steinacker 2007] 
IBeuther etal. 2010l l. Fallscheer et al. (2009) detected a mole- 
cular outflow in this source and found evidence for a large rota- 
ting structure, or toroid, perpendicular to the outflow. Evidence 
for outflow activity in this source was already found by Beuther 
et al. (2005c) and Beuther & Sridharan (2007) who found that 
there are 4.5 yum emission features at the clump edge and that 
the spectral lines of CO, CS, and SiO show broad wing emis- 
sion. The line profile of SiO(6 - 5) detected in the present study 
also indicates outflowing gas. 

The clump shows 24 fj.m emission but is dark at the Spitzer 
IRAC wavelengths (3.6, 4.5, 5.8, and 8.0 pm). Beuther et al. 
(2010) derived an SED for this source between 24 pm and 1.2 
mm, including the recent Herschel PACS and SPIRE data, and 
obtained the total luminosity of 539 (scaled to d - 3.5 kpc). 
They also found that the ratio between the total and submm lu- 
minosity (integrated longward of 400 pm) is only 11, suggesting 
that the source is very young. 

Beuflier & Sridharan (2007) detected CH3CN(6jf - 5;^) and 
CH30H(5^ - Ak) lines towards 118223 MM3, indicative of hot- 
core chemistry, and derived the abundances of 9 x 10 " and 
6.1 X 10 '° for CH3CN and CH3OH, respectively. SSK08 de- 
tected HC3N(5 - 4), only weak CH30H(7^ - 6a:) lines, and no 



CCS(43 - 32) lines towards this source. SSHIO did not detect the 
fines SO(22 - 1 1 ), OCS(8 -7), or CH30H(2i - 1 1 ) A" in their sur- 
vey; however, they detected the C2H(A^ =1-0) line. Based on 
the column densities of different species (e.g., SiO and H'^CO^), 
SSHIO suggested that 118223 MM3 is in early stage of evolu- 
tion. This conforms to the fact that the second largest value of 
^d(N2H+) in our sample (0.013) is found towai'ds 118223 MM3. 

B.6. IS0SSJ1 8364-0221 SMM1 

At the distance of ~ 2.5 kpc, the clump J 18364 SMMl is the 
nearest source in our sample. It is also the coldest ( 1 1 .4 K) clump 
of our sample. The clump is dark at 8 pm but it is associated with 
the 24-pm point source. The second largest value of 7?d(HCO^) 
in our sample, 0.012, is found towards this source. It also shows 
high cosmic-ray ionisation rate of H2 (~ 5x10""' s"'). Birkmann 
et al. (2006) studied this clump through 7 = 3-2 transitions of 
HCO^ and H'^'CO^ (beam size ~ 9'.'5). The former line showed 
blue asymmetric profile, indicating infall. They also found sig- 
nificant CO(2- 1) line wings, indicating the presence of outflows. 
The H'^C0^(3 - 2) line we have observed shows an asymmetric 
profile with a central dip and slightly stronger red peak, contrary 
to that observed by Birkmann et al. (2006) in HCO+(3 - 2). This 
difference is probably due to the larger beam size of our obser- 
vations (24"), yielding a signature of outflowing gas motions. 

More recently, J18364 SMMl was studied in detail by 
Hennemann et al. (2009). They resolved this contracting clump 
in the interferometric mm continuum into two compact cores, 
named SMMl North and South separated by 9'.'5 (0. 12 pc). Their 
positions are indicated in Fig. [1] The peak H2 column densities 
and dust temperatures were found to be 2.7 x lO^-' cm"^ and 15 
K for SMMl North and, 2.4 x 10^^ cm^^ and 22 K for SMMl 
South. Using the revised distance 2.5 kpc, the radius, mass, and 
luminosity of the northern core are 0.06 pc, 19 M0, and 26 L0, 
whereas for the southern core these values are 0.05 pc, 13 M0 
and 230 Lq. Thus, the cores within the clump have comparab- 
le sizes and masses but the southern one is associated with the 
Spitzer 24 and 70 pm sources and is more luminous. Hennemann 
et al. (2009) found that SMMl South drives an energetic mole- 
cular outflow and that the core centre is supersonically turbulent. 
On the other hand, the IR-dark core SMMl North shows lower 
levels of turbulence, but it also drives an outflow. Both the out- 
flows from SMMl North and South are rather collimated and 
their estimated ages are < 10^ yr. The HCN(1 - 0) modelling 
results by Hennemann et al. (2009) showed that the spectrum 
of SMMl South can be explained with a collapse of the core. 
They obtained an infall velocity of 0. 14 km s"' and an estimated 
mass infall rate of ~ 3.4 x 10"^ Mq yr"' (scaled to the revised 
distance). 

In summary, the clump J18364 SMMl is fragmented into 
two cores which both harbour protostellar seeds, possibly evol- 
ving into intermediate- to high-mass stars. As discussed in 
Sect. 5.4, MHD wave propagation could have played a role in 
fragmenting the parent clump. As the southern core harbours a 
24 pm source, is highly turbulent in the central region, and shows 
jet features at large distance from the driving source, it appears 
to be more evolved than the northern core. 



